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1.  INTRODUCTION 


Applications  of  the  multicomluctor  transmission-line  computer  program,  NLINE,  are  described  in  this 


manual.  Previous  reports 
NLINE  is  based. 


have  detailed  theoretical  development  and  semi-empirical  methods  on  which 


It  is  assumed  that  the  reader  knows  computer  programming  in  FORTRAN.  The  annotated  listing  of  the 
NLINE  Standard  Deck  (NSD)  given  in  section  6 represents  the  complete  NLINE  code.  The  NSD  is 
compatible  with  the  CDC  6600/6700  computers.  If  IBM  computers  are  used,  it  is  necessarj’  to  modify  the 
alpha-numeric  statements  {i.e.,  change  dimensioning  due  to  different  word  length). 

It  is  recommended  that  an  output  plotting  capability,  compatible  with  the  user’s  computer,  be  added  to 
the  NSD  for  most  efficient  output  display. 


2.  TYPES  OF  PROBLEMS  SOLVED  BY  NLINE 

The  NLINE  code  computes  currents  and  voltages  induced  by  a plane-wave  electromagnetic  pulse  (E.MP) 
impinging  on  the  conductors  and  terminations  of  a lossless  mulliconductor  transmission  line.  Computations 
are  initially  done  in  the  frequency  domain.  Time-domain  results  are  obtained  by  a numerical  Fourier 
transform  subroutine.  The  response  depends  on  the  physical  and  electrical  characteristics  of  the  transmission 
line  and  upon  the  orientation,  pulse  shape,  and  amplitude  of  the  incident  EMP.  The  transmission  line  may  be 
located  in  free  space  or  above  a ground  plane,  and  it  is  described  by  its  length,  /,  conductor  cross  section, 
propagation  velocity,  v,  terminating  impedances,  and  choice  of  reference  conductor.  Figure  I illustrates  the 
above  properties  for  a typical  line. 

The  total  number  of  conductors  m the  transmission  line  is  designated  here  by  n,  and  by  “NWIRES"  in 
the  NSD.  A convenient  reference  conductor  (e.g.,  the  ground  plane,  if  one  is  present)  is  chosen,  and 
numbered  0,  with  the  remaining  (n-1)  conductors  numbered  consecutively  as  shown  in  figure  1(a),  where 
n = 3.  The  propagation  velocity,  v,  for  the  transmission  line  is  assumed  to  be  known.  As  a matter  of 
consistency,  in  this  report  the  x-axis  is  always  chosen  parallel  to  the  longitudinal  axis  of  the  line.  The  point  of 
voltage  and  current  computation  is  thus  designated  by  a value  of  x within  the  limits  0<.r</.  Values  of 
jr  = 0 and  x = / give  results  at  the  terminations.  Figure  1(b)  illustrates  these  parameters. 

Terminations  between  conductors  are  represented  by  the  doubly  subscripted  variables,  Z,j,  in  figure  1(a). 
NLINE  is  programmed  for  all  possible  connections  between  the  conductors.  The  termination  values  are 
assumed  known,  or  must  be  measured  or  estimated  (e.g.,  assume  open,  short,  or  matched  conditions)  for  the 
line.  The  NSD  accepts  only  real  values  for  the  The  Thevenin  Equivalent  method  described  later  in  this 
report  may  be  used  with  NLINE  if  complex  or  nonlinear  loads  are  required. 

The  incident  EMP  waveshape  is  represented  in  NLINE  by  two  independent  expressions  made  up  of 
exponential  terms.  The  user  can  choose  either  when  making  an  NLINE  run.  In  addition,  the  option  of  using 
a unit  impulse  as  the  input  has  been  incorporated  into  the  code.  The  first  of  the  exponential  pulse  forms  is 
designated  as  Pulse(l).  This  pulse  is  the  familiar  two-exponential  form  given  below. 

PulsefI).  Two-exponential  representation 


E\(t)  = A (e-''  - e-^')  / M 


where  £'(0  is  the  y-component  of  the  incident  EMP  electric  field  in  volts/meter. 

/q  - In  (a2^®i)  ^ («2“  ®i)  zero-to-peak  rise  time  of  the  pulse. 

' Sidney  Frankcl.  Cable  and  Multiconductor  Transmission  Line  Analysis,  Harry  Diamond  Laboratones  TR-09M  (28  June  1974). 

^ Jams  Klcbcrs.  Electromagnetic  Pulse  Coupling  with  l.ossless  Multiconductor  Transmission  Lines.  Harry  Diamond  Laboratories  TR-171 1 (June 
1975). 
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,4  — peak  amplitude  of  the  pulse  in  volts  meter. 

M = normalization  factor  which  makes  the  pulse  reach  amplitude,  A.  for  all  a,  and  Oi. 
a I and  uj  are  pulse-shaping  constants. 

The  second  ol  the  e.xponential  pulse  forms  is  designated  as  Pulse(  I ).  This  pulse  is  made  up  of  several 
exponential  terms,  has  a zero  time  derivative  at  time  equal  to  zero,  and  can  be  made  to  have  a zero  cross- 
over at  designated  late  times.  Pulset  I)  is  useful  in  representing  EMP  simulator  fields  such  as  those  generated 
by  AESOP  (Army  Electromagnetic  Pulse  Simulator  Operations)  and  TEMPS  (Transportable  Electromagnetic 
Pulse  Simulator)  simulators  at  HDL.-* 

Pulse(-!),  llme-donuiin  expretcsion  for  the  electric  field* 

£'(/)  = (u-/>)(l  -I-  ^ + {b  + d){\  + hl)e''" 


- u(l  + ml)e  — (J(\  + pt)e  ^ 

where  the  pulse-shaping  constants  are  given  below  with  their  corresponding  FORTRAN  variable  equivalents 
in  NLINE. 

m = Bl,  A = B2,  f!  = B3,  p = B4,  a = B5,  b = B6.  d = B7. 

The  basic  waveshapes  of  Pulsed ) and  Pulse(  I)  are  shown  in  figure  2. 


PULSE(I) 


Figure  2.  Basic  pulse  shapes  of  the  input  EMP  to  NLINE. 


% 


3£.  Patrick.  S.  SooHoo.  Transportable  Electromagnetic  Pulse  Simulator  (TEMPS)  Preliminary  Field  Mapping  Report.  Ham  Diamond 
Laboratories  TM-74- 1 5 (October  1974). 

* After  a pulse  representation  bv  R.  F.  Gray.  HDL. 
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3.  OUTLINE  OF  EQUAIIONS  ISED  IN  NI.INE 


NLINE  IS  based  on  a solution  for  the  I I.V1  mode  currents  and  voltages  induced  by  the  distributed 
external  electromagnetic  fields  along  the  conductors.  The  equations  governing  the  values  of  the  current.  /, 
and  the  voltage  y.  along  the  axial  x-direction  of  the  transmission  line  are  given  below. 


clj 

Jx 


^ + ;<a  r K = //' 


(3) 


jy 


+ yw  LL^  E! 


dx  - (4) 

where  the  underlined  quantities  denote  matrix  form.  The  matrices  are  of  the  order  (n-1).  Solutions  for  the 
voltages  y_  and  the  currents  / are  derived  in  Frankel,'  and  are  given  in  the  form 


/ = l'(^'costf  + yl  sin»)s‘'^(/)  + ]£(.v 


0 


r = (icos^  + ye'sin#)5  'K(!)  + U_(x)  . 


(5) 


(6) 


The  FORTRAN  version  of  equations  (5)  and  (6)  is  in  the  frequency  domain  and  combined  with  the 
numerical  Fourier  transform,  constitutes  the  NLINE  computer  code. 


£ 

L 


The 


the  unit  matrix  order  of  (n-1). 

number  of  conductors  in  the  transmission  line. 

kx 

u / V 

radian  frequency  in  rad  s. 
transmission-line  propagation  velocity. 

Maxwell’s  capacitance  coefficient  matrix  for  the  conductor  cross  section. 

inductive  coefficient  matrix  obtained  from  the  C matrix,  unit  matrix  I_.  and  the  line  propagation 
velocity  v by  the  following  relation 


= -yi 


driving  functions  in  equations  (3)  and  (4)  are  given  by 

/£'  = (wi.  H, 

r = (fi.  El E„^,y . 


0) 


(8) 


' Sidney  Frankel.  CabU  and  Muluconductor  Tranimission  Ltne  Analysis.  Harry  Diamond  Laboratoncs  TR-091-1  (2S  June  1974). 


..^ 


where  >'  denotes  transpose,  converting  tC  and  to  column  matrices. 
The  /'*'  elements  of  equations  (2),  i = 1,2 (n-  I i are  defined  by 


//'  = juC‘l:  y{u)e 


jiwl  - k..x  casv) 


E,  = juL,H,{u)e‘  . 

Terms  in  equations  (9)  are  outlined  below. 
i “ time  in  seconds. 

- angle  of  EMP  incidence  shown  in  figure  1(b). 

Ey(u)  “ the  resultant  y-component  of  the  external  electric  field,  in  volts  meter. 

H y(u)  = the  resultant  z-component  of  the  external  magnetic  field,  in  amperes/ meter. 


The  magnetic-field  coupling  parameter  is  defined  in  transpose  form  as  follows: 

4'  = (T^  L\ T'.,)V 

ij  is  an  input  parameter  to  NLINE,  and  is  defined  by  Frankel*  and  Klebers.^ 


c'  = in  c\ 


r'  is  computed  in  NLINE  from 


I n-  I 

r'  = - — y L'c 

^ I- > 


where  is  the  permeability  of  free  space. 


The  electric  and  magnetic  fields  are  integrated  over  the  parameter  .t  in  the  following  relations. 

LLix)  = J cos|A:(x-01  ~ sin[*(xr  - 01  ) 

o 

IK(X)  = T/ cos(A:(x--01  - ;£'(0  sin|/c(a:-01  } 


Z = V L 


Y = V C . 


’ Sidney  Frankcl.  Cable  and  Muliiconducior  Transmission  Ijne  Analysis.  Harr\  Diamond  I,aboratoncs  TR-09I  - 1 (28  June  1974). 

* Jams  Klebcrs.  Eleciromaj(nelic  Puhe  Coupling  ^ith  l.ossless  Muliiconducior  Transmission  Lines.  Harrv  Diamond  l^biiralones  TR  - 171 1 (June 
1975). 


\ 


k 

k 


h- 


Other  terms  in  equations  (3)  and  (4)  are 

S=(^p'+  t')  coskl  + >(i  + t'E!)  sin  kf 


(16) 


t'  = ZLT 


(17) 


p‘  = m' 


(18) 


K(h  = z[]vu)  - ru[f)] . 


(19) 


The  load  impedances  ot  the  transmission  line  are  introduced  in  the  equations  through  admittance 
matrices  and  Y^,  where  the  superscript  “i"  refers  to  loads  at  x = 0.  and  the  superscript  "o"  refers  to 
loads  at  ,e  = /. 

Elements  of  K'  and  are  defined  by 


1 

! 

1 


>':i  = 


n - 1 


(20) 


= 


n ■"  1 

2 = / 

/-I 


(21) 


where  )^J"  is  the  load  admittance  between  the  /'*'  and  j'^  conductors  of  the  transmission  line  at  the  .v  = 0 
end  in  equation  (20),  and  at  the  x = f end  in  equation  (21),  respectively. 

In  each  case,  the  impedance  between  the  i'*'  and  conductors  is 


The  Zj  values  are  input  parameters  to  NLINE. 


(22) 


4.  DEFINITION  OF  NLINE  INPUT  PARAMETERS 


The  input  parameters  are  defined  and  listed  in  sequence,  as  they  are  read  from  data  cards  in  NLINE. 
The  two-exponential  Pul.se(l)  is  read  by  data  cards.  The  parameters  for  Pulse(-l)  are  changed  within  the 
NLINE  deck.  Setting  the  FORTRAN  variable  IMP  = I in  the  NSD  implies  use  of  Pulse(l)  for  the 
computation.  Setting  IMP  = - I implies  use  of  Pul.se(  I)  for  the  computation. 

The  FORTRAN  variables  read  by  data  cards  into  NLINE  in  the  following  order: 

A - amplitude  of  the  incident  EM P in  volts/meter. 

A1  - a I in  equation  (I). 

A2  « a 2 in  equation  ( 1 ). 


j 


1 

I 

I 


! 

1 


10 


The  FORTRAN  statements  as  thc\  appear  m NUNT-.  tor  the  above  and  subsequent  variables  are; 


READ(  5,5000  ) A.  A1  . A2 

5000  FORMAT (3E  10  0) 

PFIRCF-N  = percent  of  tree  space  propagation  velocitv  observed  in  the  transmission  line. 

READ(5 . 5001 )PERCEN 

5001  F0RMAT(F  5 0) 


COMNT  - user's  comment  (up  to  80  characters) 

READ(S . 505) (COMNT (MN) , MN=1 . 8 ) 

505  FORMAT(8A10) 


MSETS  = number  of  data  sets  where  values  for  the  transmission-line  terminations,  and  or  angle  of 
incidence,  and  or  line  length  ma\  be  parametricallv  changed,  while  holding  reference  conductor,  cable  cross 
section,  and  a.xial  orientation  with  respect  to  the  incident  LMP  constant. 

NWIRES  = total  number  of  conductors  in  the  transmission  line.  NLINF.  is  dimensioned  so  that  the 
limit  for  NWIRES  is  I < NWIRES<  I 1 . 

READ( 5 , 500 )MSETS . NWIRES 
500  FORMAT(2I5) 

C.MAXWL  = Ma.xwelTs  capacitance  coelTicientsT  for  the  transmission-line  conductor  configuration.  .A 
reduced  matrix  is  used  because  the  off-diagonal  terms  are  symmetrical.  These  coefficients  can  be  measured  or 
computed  by  various  methods.  A program,  PfNTL  (Capacitance  and  Magnetic  Field  Coupling  Code),  is 
included  in  appendix  .A  for  one  computation  method. 

The  capacitance  matrix  is 


(CMAXWL) 


‘C„  C,2  . 

• ■ c,„- 

C22 

..  C2„ 

Svmmetric 
terms  ni't 

punched  or 
data  cards. 

c;.! 

For  example  n = 3 for  a three-conductor  transmission  line.  The  CMAXWL  array  would  be  punched  in 
the  order  given  below: 


C|2.  C|,.  C22. 


where  the  symmetrical  off-diagonal  terms  C21.  Cj,.  and  C',2  are  not  required.  In  general,  for  n conductors. 
n{n-¥  \ )/l  terms  are  read  into  the  CMAXWL  array. 

NCAPS=NWIRES*{NWIRES+1 ) /2 
READ(5 , 501 )(CMAXWL( J) . J+1 .NCAPS) 

501  FORMAT(8E10  0) 

It  is  important  to  note  that  the  CMAXWL  array  is  converted  to  the  full  rt  by  n-capacitance  array  in 
NLINE.  However,  the  1 row  and  1 column  are  deleted  from  the  array,  after  the  i conductor  has  been 
designated  as  the  reference  conductor. 

^S.  Ramo  and  J.  R.  Whmncrv.  Fields  anii  H aves  in  Minicrn  Radio,  .lohn  Wiley  & Sons.  Inc..  New  \'ork  ( 1%4),  264-265. 


'-}i^ 


1 


\ 


For  example:  n ^ 3.  the  Maxwell's  capacitance  arras  generated  from  CMAXWI.  m NLINK  i.s: 


C„ 

(-12 

C|3 

c = 

^^21 

C22 

(^2.1 

-('31 

Q2 

(23) 


For  example,  if  conductor  number  three  is  chosen  as  the  reference  (i  = 3),  it  is  then  redesignated  as  the 
0'*’  conductor,  and  the  reduced  capacitance  matrix  used  for  the  compulations  becomes 


C'  = 


'c„  c,r 

LCii  C22— 


(24) 


Hence,  the  values  for  C,j,  C,,.  and  Cj,  are  not  needed  if  the  choice  of  reference  conductor  remains  fixed 
for  all  computations. 

ELI  = the  magnetic-field  coupling  parameter  array,  L‘,  The  (NWIRES-1)  terms  are  written  into 
NLINE. 

READ(5,503)(ELI(J)  , J=1  , (NVKIRES-1  ) 

503  FORMAT(8£10  0) 

The  magnetic-field  coupling  parameters  can  be  measured,  or  computed.''^  The  PTN'l  L program  listed  in 
appendix  A performs  this  computation. 

IQ  = the  integer  number  designating  the  chosen  reference  conductor. 

READ(5.504)IO 

504  FORMAT(I5) 

Z1  = Z*",  the  array  of  terminating  impedances  at  the  ar  = 0 end  of  the  line. 

ZO  = Z*",  the  array  of  terminating  impedances  at  the  x=!  end  of  the  line. 


After  IQ  is  chosen,  the  reference  conductor  is  renumbered  as  zero  (0),  and  the  remaining  conductors  are 
renumbered  consecutively  1 through  (n-I).  The  like  subscripts,  lY,  designate  impedance  connected  between 
the  i conductor  and  the  reference  conductor  (0),  The  subscripts,  ij.  designate  impedances  between  the  / 
and  y conductors.  The  values  for  ZI  and  ZO  are  read  row  by  row,  left  to  right.  For  example:  n = 3,  the 
sequence  of  zjj  stored  in  ZI  is 


Z,2.  ZzV  zi) 


N2-(NWIRES-1  )*»2 

READ(5 , 501 )(ZI(J) , J=1 . N2 ) 

READ(5 , 501 )(ZO( J ) , J=1 ,N2 ) 

501  FORMAT(8E10  0) 

PSI  * ly,  the  angle  of  incidence,  in  degrees,  shown  in  figure  1(b). 

TL  » /,  line  length,  in  meters. 

REAO(5.502)PSl,TL 

502  FORMAT(2F10  0) 

* Sidney  Frankel.  Cable  and  Muluamducior  Transmission  Line  Analysis.  Harry  [>iarm>nd  l-aht^ra lories  TR-Wl-l  (28  June  1974). 

* Jams  Klebers.  Eleclronuigneiu  Pulse  Coupling  wiih  Lossless  Muhu'onduaor  Transmission  Lines,  Harrv  [>iamond  1-ahoraloncs  TR  - 171 1 (June 
1975). 
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The  remaining  paramelers  under  user  conirol  in  NLINE  are  punched  and  changed  within  the  deck  itself. 
These  parameters  are  listed  below. 

IMP  = the  parameter  that  determines  which  input  HMP  field  is  used.  IMP  = 0 calls  the  unit  impulse, 
1MP“  1,  calls  Pulsed),  and  1MP  = -I  calls  Pulsei  1), 

Bl,  B2,  B3,  B4,  B5,  B6,  and  B7  are  the  Pulset  1)  parameters  listed  in  equation  (2). 

IDATE  = any  integer  number  set  bv  user  as  an  identification  number  for  the  computation. 

NXTl.ME  = number  of  desired  axial  .x-positions  at  which  the  computation  is  made. 

X87  = fraction  of  the  line  length  at  which  the  NXTIME  computations  are  to  be  made.  For  e.xample,  if 
computations  are  desired  at  ,v  = 0.  .v  = 0.5/,  and  x = I,  holding  all  other  line  parameters  constant,  the  user 
sets  NXTIME  = 3 and  X87(l)  = 0.  , X87(2)  = 0.5,  X87(3)=1.0. 

NCOUNT  = number  of  frequency  points  at  which  the  computations  are  made. 

NTIME  = number  of  time  points  at  which  the  computations  are  made. 

These  parameters  can  be  changed  to  adjust  NLIN’E  run  time  and  Fourier  transform  accuracy.  Maximum 
allowable  values  for  NCOUNT  and  NTIME  are  600. 

5.  ABBREVIATED  FLOW  DIAGRA.M  FOR  THE  .MAI.N  DECK 

A flow  diagram  for  the  NLINE  main  program  is  listed  in  this  section.  Only  the  basic  sequential  events 
are  included. 
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J 


J877— 1, 

. NXTIME 


FREQUENCY- 

INDEPENDENT 

COMPUTATIONS 


f BEGIN  ' 
FREQUENCY 
V DOMAIN  J 


NCOUNT 


USER  SPECIFIES  NUMBER  OF  FREQUENCY 
POINTS  AT  WHICH  COMPUTATION  IS  MADE 


begin  FREQUENCY- 
DOMAIN 
COMPUTATIONS 


CURRENT  AND  VOLTAGE 
MATRIX  COMPUTATIONS 
COMPLETED 


BEGIN  TIME- 
DOMAIN 

COMPUTATIONS 


USER  SPECIFIES  NUMBER  OF  TIME  POINTS 
AT  WHICH  COMPUTATION  IS  MADE 


FOURIER  TRANSFORM 
FROM  FREQUENCY  TO 
TIME  DOMAIN 


PRINT  TIME- 
DOMAIN 
RESULTS 


OPTIONAL:  USER  MAY  PLOT  TIME  AND 
FREQUENCY  RESULTS 


6.  NLINE  STANDARD  DIXK  (NSD)  LISTING 


This  section  includes  the  listing  for  the  entire  NLINE  computer  program,  including  the  mam  program 
and  all  subroutines  and  functions  required  for  an  NLINE  run.  A considerable  amount  of  comment  cards  has 
been  added  for  the  user’s  benefit.  Three  types  of  comment  cards  have  been  placed  in  the  deck.  A plain  "C” 
in  the  listing  left  margin  designates  location  of  an  informative  comment  for  the  user.  A “CR”  in  the  margin 
locates  a position  in  the  deck  where  data  are  read  in  from  data  cards.  Finally,  a “CVJ"  in  the  margin  locates 
a point  in  the  deck  where  the  user  can  modify  a particular  parameter  within  the  NLINE  deck.  ITiis  listing 
was  tested  on  the  CDC  6700  computer. 
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PRCOKAV  MlNt  (INPUT, OUT  PJT,T& PE  5»INPJT,JUK,TAPE2  = A4K, 
lTflPE6=CUTPUT ,TAPt 1 ) 

C 

C 

C MINE  M^LTUDNOUCUR  TkftNS'^ISSUN  LINE  COKPUTfR  CCCE.  JUNE  1976. 

C THIS  wr^K  hAS  SPONSORED  3Y  THE  DEFEN5E  NUCLEAR  AGENCY  UNDER 

C SJETASK  R99.AXEPD7S. 

C DEPART«rNT  OF  THE  ARMY 

C HARRY  OIAMCNC  LASCRATCRIES 

c 2&C0  plbDer  mill  RD 

C ADELPHI  , MD  2C7E3 

C 

C 

DIMENSION  CXMC(60C,1D),CXAT(SOO,10),YXMG(6UO,10),VXAT(600,10) 
DIMENSION  E2TTC6QOI,CSrtC(6OOI,EERK(60C).FKEK(bCO),TlME (60C I , 
lCTIM(600),VTIh(63C),H2FTM(6O0) 

DIMENSION  CMAXWLi  lUOl  ,71  { 100)  ,2C(  IOC)  ,Y0(  IOC)  ,YI  ( 1001  ,P0UC0)  , 
lPMlOO),YnOO),C(lCO),Z(100),PP(lOO),SOnOO),C(lCC),CAP{100), 
2XPtl  00) 

0 IMfc,\iSI2N  X37(10)  ,£  LI  (ID)  ,VXMA&(  10)  rCXMAGdO) 

DIMENSION  CuMNT  13  ) 

Complex  e2TEssct»oci  ,temd(6j0) 

COMPLEX  110(600) 

COMPLEX  RC(lCO),CS(ljO),Y3C(lOO),CnMP2(100),PVX(10C),PlX(100), 
IVX(IOD) ,CX(100) 

CCPPLEX  JX(lU),WX(10)iUL(lD)thL(10),XKL(10) 

COMPLEX  2UZ9,SJ,H7FREE 

C0MMJN/XYr/KTCNT,NTH,TL,»51tX 

CC.M'»0N/K&L/V,lMP 

commjn/xsav/cap 

C:>MM3N/FPAR/A  1 ,A?  ,A 

CCKMUN/YCDLL/Pl,32,b3,BA,B5f66,P7 

data  P I E .UM  »kAD/3  .1^1 S92  7 , 1 .2S663706E -6, .C 1 7A532S/ 

CALL  SYSTEMCOC, 0,1  ,1,1, -1,-1  ) 
call  SYSTEMCni^.C.l,l,l,-l.-l) 

C 

C Two  EXPQNENTUL  PULSE  PARAMETERS  ARE  READ  HERE. 

C REFER  TO  COMPLEX  FUNCTION  HZFREE. 

CR 

REA0(5,S000)A,A1,A2 

5000  FCRMATOEIO.O) 

C 

C PERCeN*PERCENT  FREE  SPACE  VELOCITY  OF  LIKE 

CR 

REA0(S,5001)PERCEN 

5001  FDRMAT(F5.0) 

V*PERCEN»3.E8 
SJ*CMPLX(0.,1  .) 

C 

CU  USER  SETS  DESIRED  VALUE  OF  (IMP)  HERE. 

C IMP  SPECIFIES  TYPE  OF  INCIDENT  PLANE  WAVE  PULSE  SHAPE. 

C IMP«0,  HZFREE^UNIT  IMPULSE 

C IMP*-l,HZfREE»SPECIAL  A-EXPONCMT I AL  PULSE 

C ,H2FREE*TW0  EXPONENTIAL  PULSE 

C three  CAROS  ARE  INCLUDED  HERE  FOR  THE  THREE  VALUES  OF  IMP. 

C PLACE  DESIRED  VALUE  OF  IMP  LAST  IN  CARO  SUCCESSION 

C 

1 MP*0 
IMP«1 
IMP*-1 
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r 


lF(lliP.6Q.-llGa  TO  110b 
GO  TO  1107 
1106  CONTINUE 
C 

C SAMPLE  VALUES  FOR  THE  4'E XPQNENT 1 AL  PULSE  PARAMETERS  ARE  GIVEN  HERE. 

C REFER  TO  COMPLEX  FUNCTION  HZFREE. 

C 

A«1  .E2 


I 


k 


k 


110b 

1107 


C 

cu 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

CR 


SOS 

602 

C 

C 

C 

c 

c 

c 

c 

CR 


bOO 


c 

C 

c 

CR 


bOl 

c 

c 

C 

CR 


503 


B1«7.E6 

62«6.5E6 

63«0. 

B<»*0. 

B5*b.Eb 

B6=b.ER 

B7=0. 

T=0. 

DTAU=b.E“lO 
AMPT0*0  . 

DC  110b  NJ03«I,300 
T*T*OTAU 

CALL  TEMPLT(T*ETST,A,Al,A2,IMPI 

IF  (ETST.CT.AMPTfi)AMPTQsETST 

CONTINUE 

A=A*A/AMPTQ 

CONTINUE 

IRUNO^O 

USER  SETS  DESIRED  VALUE  OF  RUN  IDENTIFICATION  NUMBER  HERE. 
lOATE^RUN  IDENTIFICATION  NUMBER. 

FOR  EXAMPLE.  DESIGNATE  DATE  OF  RUN  AS  FOLLOWS. 

IDATE=04066 

IDATE^IIJJK 

I I*DAY, JJ=MDNTH.K»YEAR 

I0NTFN*I0ATE*1000 

IDNTFN=I0NTFN4IRJNQ 

IDNTFN^l.D.  NUMBER^IIJJKLLL.LLL^RUN  NO. 


USER  GENERATED  COMMENT  OF  FORTRAN  CHARACTERS  READ  IN  FROM 
COLUMNS  1 through  80  ON  DATA  CARD. 

REA0(5,505  MCOMNT  (MNl  .MN«  1.61 
FORMAT(BAIO) 

WRlTE(6»b02MC0MNTCMN).MN»l«e) 

FORMAT!  1H1.8A10//  ) 


MSETS’NUMBER  OF  DATA  SETS  WHERE  VALUES  FOR  THE  TRANSMISSION  LINE 

TERMINATIONS  UI.ZOI.  ANO/OR  ANGLE  OF  EM  NAVE  INCIDENCE  CPSIl. 
ANO/DR  TRANSMISSION  LINE  LENGTH  (TL)  MAY  BE  PARAMETRICALLY 
CHANGED.  WHILE  HOLDING  REFERENCE  CONDUCTOR.  CABLE  CROSS  SECTION. 
AND  AXIAL  ORIENTATION  WITH  RESPECT  TO  INCIDENT  EM  FIELD  CONSTANT. 


Nk IRES=NUM6ER  OF  CABLE  CONDUCTORS  INCLUDING  REFERENCE  06 

REAO(b.bOO)MSET5.NHlRES  Ob 

F0RMAT(21bI  06 

NCAPS=NWlRES«(NNlRES«ll/2  06 

NP^NhIRES-1  06 

N2=NP*NP  06 

MAXWELL  CAPACITANCE  COEFFICIENTS  FOR  THE 
TRANSMISSION  LINE  ARE  READ  AT  THIS  POINT. 

READ(b,b01)(CHAXwLU)»J>l.NCAPS)  06 

FCRMATtSElO.C)  06 

MAGNETIC  FIELD  COUPLING  PARAMETERS  FOR  THE 
TRANSMISSION  LINE  ARE  READ  AT  THIS  POINT. 

READ(b.ba3MELIIJI.J»l.NP)  06 

FCRMAT(8E10.0)  06 


I 


1 


i 


:? 
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C tC^RtFfcittNCE  CONDUCTOR  NUHbtR.  ONE  3F  TmE  INnIRES)  NUHBER  OF  CONDUCTORS 

C IN  the  TRANSNISSION  UNE  IS  ChOSEN  AS  THE  REFERENCE- 

C« 

REACIS.SD^MQ 

SON  fCRNATUlSI  0^ 

00  S6  JJS8>1.HSETS 
C 

c READ  TERMINATION  RESISTANCE  VALUES  AT  X-0  ENO  (Zll  AND  l«TL  END  IZOI 

C OF  THE  TRANSMISSION  LINE.  TL«LENCTH  OF  THE  TRANSMISSION  LINE. 

CR 

REA0(S,S01MZ4(J)  ,J«1,N?) 

READIS.SOl  MZOI  J)  tJ^ltN?) 

CALL  RAR(Z1.2HZI.NP> 

CALL  RAR(Z0,2HZ0*NP  I 
C 

C ANGLE  OF  EMP  INCIDENCE  (PSII,  AND  TRANSMISSION  LINE  LENGTH  ITLI 

C ARE  READ  HERE. 

CR 

READ(S,502)PSI.TL 
S02  F0RMAT(2F10.01 

CALL  RVEC (ELI .2HLE .NP I 
C 

CU  USER  SPECIFIES  (NXTIME)  AND  XBTdl.AS  DEFINED  6EL0M*  HERE. 

C NXTIME«NUMBER  OF  DESIRED  AXIAL  POSITIONS  OF  COMPUTATION 

C ALONG  THE  CABLE. 

C X87(I ) ^FRACTION  OF  TOTAL  LENGTH  OF  CABLE  AT  MHICH 

C COMPUTATION  IS  MADE#  1 «1 • 2 f • . ..NXT IME . 

C COMPUTATION  IS  MADE  AT  THE  AXIAL  POSITION*  X«Xe7(l)*TL*  MHERE 

C TL*  LENGTH  OF  THE  CABLE. 

C 

NXTIME^l 
X67U  1-0. 

X87(2)*0.2S 

C 

DO  66  J877«l .NXTIME 
X*X87( J877»*TL 
CALL  rOYI (Y0,Y1 ,ZO.II #NP I 
CALL  RAR(Y0t2HY0»NP ) 

CALL  RAR(Y1*2HY1 .NP I 

CALL  POPI  (NrIRES,  1Q,Y0*YUPQ,P1,Y,Z,CMAXNL) 

CALL  RARiP0*2HP0,NP ) 

CALL  ftAttfPl 

C 

C*»**8EG1N  FREQUENCY  DEPENDENT  COMPUTATIONS. 

C 

CU  USER  SPECIFIES  INCOUNTI*  AS  DEFINED  BELOW*  HERE. 

C NCOUNT»NUMBER  OF  FREQUENCY  POINTS. 

C 

NCQUNT>bOO 

C 

C OMEGl  FILLS  THE  ARRAY  FREK  WITH  THE  RADIAN  FREQUENCY  VALUES. 

C SEE  SUBROUTINE  OMEGl  FOR  FREQUENCY  GENERATION  SCHEME. 

CALL  0MEG1IFREK*NC0UNT) 

00  87  LL8*1,NC0UNT 
W«FREK(LL8) 

BETA»W/V 

theta=beta»tl 

cc*cosithetaj 

SN*5IN(THETA) 

CALL  DADO  (P0*P1 *PP,N2) 

CALL  SMLT  IPP *CC *NP,PP I 
CALL  SMLT(P0*SN*NP*S01 
DO  56  JXXsl*N2 
56  XP( JXX|=PI(JXXI 

CALL  MPRD(Sa*XP,Q*NP*NP*0*0*NP) 

CALL  XIDENT(NP.SO) 

DC  10  J«1 *N2 
10  RC(J)«CMPLX(SO{J),0.) 

CALL  SMLT(SO,SN*NP,SOI 
CALL  DADD  (S0*0f0*N2 I 
00  13  J>1*N2 
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c 

C CS*5  MATRIX 

C 

13  CS(  J)<CRPLX(PP(  J)  ,C(JM 

CALL  COELC  (RC  iCS.NP  ,NP,1  .f  >10  ,UR  ) 

C 

C RC«S  INVERSE  MATRIX 

C 

CALL  UwOPMX,h,2.PSl,UX.MX,y,NP,fin 
CALL  UWOFX4TL  ,l  .PSl  ,JL,ML  ,Y  .NP,ELn 
DD  11  J«l  ,N2 
CO«PZ( Jl»CMPLX(Z{ J»  ,0.1 
11  TOC(J)»CNPLXi-YOU)  ,0.1 

CALL  CHPRD(TOC,Ul ,C S,HP «NP , 0 • 0,1 > 

CALL  CMAODCML ,CS,CS,NP) 

CALL  CMPtO(CORPZ,CS  .IKL .MPtHP  ,0,0, 11 
C 

C XKL*K(U  VECTOR 

C 

CALI  MN(X,6ETA,PVX.PIX,Y  ,P(,NPi 
CALL  CHPROERC ,XKL ,C  S,NP, NP, 0,0,1) 

DO  IS  K-l,NP 
IS  CX(K)-CS(K1 

CALL  CN'PR0lPVX,CS,XKL,NP,NP,0,0,l ) 

CALL  CHADD(XKL,UX.VX,NP) 

2U2ft*C£  XP lSj*e€TA*COS(PS I •RAO )*TL ) 

CALL  CSHLT(VX,ZUZ6,NP,VX) 

c 

C VX*V0LTA&£  vector  at  X 

c 

CALL  CHPR0(P1X,CX,XKL,NP,NP ,0,0,1) 

CALL  CMAOO  (XKL  ,WX  ,CX,NP) 

CALL  CSMLT(CX,ZUZ8,NP,CX) 

C 

C CX*CURRENT  VECTOR  AT  X 

C 

DO  1<»  J*1,NP 

CALL  SMAGICXT J) ,CXN&(LL8, J) ,C  XAT(LL8, J)1 
lA  CALL  SMAG(VX(J)»VXN&(LL8,J),VXAT(LL8,J)) 
87  CONTINUE 
C 

C»»A*fMD  FREQUENCT  DEPENDENT  COMPUTATIONS. 


C FREDUENCY  OOHAIN  RESULTS  ARE  PRfNTfO  OUT  8r  SUBROUrZNE  (HEAD)  AT  fVEtr 

C J'Th  FREQUENCY  VALUE.  J*1  PRINTS  OUT  ENTIRE  FREQUENCY  ARRAY  OF  RESULTS. 

CU  USER  SPECIFIES  VALUE  OF  (J)  HERE. 

C 

J*30 

CALL  HEA0(0,J,FREK,CXM&,CXAT, VXHG,VXAT,C3MNT,TL,PSI flQ, 
1NWIRES,10NTFN,NP,NC0UNT ,X) 

C 

C«***BEGIN  TINE  OOHAIN  CUHPUT AT  I QNS  . 

C TIME  ARRAY  IS  GENERATED  HERE. 

C the  uses  may  NISH  TO  MODIFY  FREQUENCY  AND  TINE 

C ARRAYS  TO  ACHIEVE  HOST  ACCURATE  FOURIER  TRANSFORM. 

CU  USER  SPECIFIES  (NTINE),  AS  DEFINED  PELON,  HERE. 

C NTINE-NUN8ER  CF  TINE  POINTS. 

C 

NTINE=600 

TAU=TL/V 

DDTAU'IOO. 

DTAU=TAU/DOTAU 
TIMCni^l  .E-12 
I TMtt  = l 

DC  1003  1TIM*2,S0 
ITIMd=inM8*l 

1003  TI»*E(ITIM)=TIME(ITIH-1)*.1»0TAU 
NTMME'NTIME-1TIH9 
ITIN8»iTINfi4l 
DC  713  IT INsiTIMB  ,NTNN| 


713  T|mm  n IM  I*  UHt  II  T IN-n  ♦OUU 
MlMt 

DC.  8-9  NTHsl.KP 
OU  90  JJ«1  .SCOUNT 

90  TfpPUj)«VX*10(JJ,NTH  »®:tXP  (S  J*VXAT(JJ,NTM  1/57.2958) 

C 

C FCURICR  transform  FROH  FREQUENCY  TO  TIME  DONAIK  PERFORMED  HERE. 

C USER  **Ar  hI5H  TO  SUBSTITUTE  PERSONAL  PREFERRED  SUBROUTINE  HERE. 

C SUBROUTINE  TRANS  IS  USED  HERE  FOR  THE  FOURIER  TRANSFORN.  j 

C i 

CALL  r RANSiTEMP.  vr  I N,FRE  K .nNE.NCOUNT  .NTIME,!  .,0)  i 

DO  91  jj>  1 .NCOUNT 

91  TENPUJ)>CXMC(  JJ,  NTH)  •CEXP(SJ*CXAT(JJ. nth  1/57.2958) 

C 

C SUbROUTINE  TRANS  IS  USED  HERE  FQR  THE  FOURIER  TRANSFORN. 

C 

CALL  transitenp.ctin.frek.tine.ncount ,NT1ME,1 .,0) 

c 

C TINE  OLNAIN  RESULTS  ARE  PRINTED  OUT  8Y  SUBROUTINE  (HTIMI  AT  EVERY  J-TH 

C TINE  VALUE.  J>1  PRINTS  OUT  ENTIRE  TINE  ARRAY  OF  RESULTS. 

CU  USER  SPECIFIES  VALUE  OF  ( J)  HERE. 

C 

J*20 

call  HTTN{CDNNT»lDNTFNtVTlN«CTlN,TlNE«NTINE»J) 

C 

CU  USER  NAY  PLOT  FREQUENCY  AND  TINE  DDNAIN  RESULTS  FOR  RESPONSE  AT  EACH 
C CONDUCTOR  HERE.  USER  NUST  WRITE  OhN  PLDTTIHC  INSTRUCTIONS  CONPATIBLE 

C WITH  LOCAL  CONPUTFR  . 

C 

89  CONTINUE 
C 

C**«*  END  tine  DDNAIN  C ON  PUT A T 1 ONS  . 

C 

86  continue 

DO  2221  1 LXJK  «1 .NCOUNT 

22211  E 2TE  S5  UXUK )»377.«HZFREE (FREKCLXUK ) .A  *A1 . A2. INP) 

CALL  TR4NS(f2Tf  SS.E2rT,  FRF»l,  TINE  .NCOUNT.N  TINE  ,1.,0) 

C 

C TINE  DCNAIN  ELECTRIC  FIELD  USED  AS  THE 

C DRIVING  SOURCE  IS  PRINTED  AT  THIS  POINT. 

C 

DO  22212  J>1  ,NT  HE 

22212  CALL  T E HP L T ( T IN E U ) tHZF TN | J ) , A, A 1 • A2 , INP ) 

MR  ITE  (t  .700)  IDNTFN,  (CONNT  HN)  »NN>1  ,e)  flDNTFN 

700  FCRHAT  ( 1M1.5HI  .D.  *,  l9,nX8A10  rlCXSHI.O.*,  19//) 

HRITE(b*7Cl)  J 

701  FORHATUH  ,85H1NC1DENT  ELECTRIC  FIELD  FRQN  EXACT  EQUATION  AND  FRON  { 

1 FOURIFR  TRANSFORM  IN  VOL TS/N f T E R / ) 

HR|TE(6t702l 

702  FORNATdH  »AX9HT  1 NE  ISCC  1 • 7ISN8IACT  .MtMTRAASPOM) 

WR1TE(6,703MTIIIC  UI»HZFTNU)  »E2TT(JI  ,J«l  .IITIflE  ,20) 

703  F0RNAT(1P3E1A.A) 

STOP 

ENO 


SU8«0JT]N£  3<<E&1(M(k,NC30IITI 

c This  &(/6«aur{Aif  the  4**4r. 

omeNsfaN  MfKibaoi 

P|(t3.l<>lS927 

C HPO^NJ'^BrR  3F  FREQUENCY  POINTS  PEN  3EC«0E. 

NPO«  l?0 

NOC*NCnjNT/NPD 

C 0H3HEST  FREQJENCV  VRIJE  . 

0-1.7^ 

0»1.E7 
0*1. F<* 

FR«3 

DN*D/FL0*r(NP0l»10. 

u o 

03  i J*l  ,NOC 
03  ) L*l tNP3 

FR«FR»C« 

J FRF^fUl  .F» 
i DP«10.*DN 
«E  rjRN 
ENO 
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C Vf  ( I V(C  .lABl  ,N(*I 

CVE:  PMMS  3jI  iHt  VECTOR  IVKI  H««E0  URBil 

CORPIEI  VEC( i > 

MR  IT  E Ifc ,9 ILABl 
9 F3*«RTnH  ,lkil 
hR1TE(6.  lOH  VEC  ( J) , 

10  F0RR*T(10ni.3l 
RETJRN 
ESO 


Subroutine  KiofNTiN.Ru 

XIOENT  GENERATES  AN  N ftT  N UNIT  MATR]l. 
CALLED  FRON  F IN.^ORI 
D1 NE  NS I3N  R 1 ( 1 » 
l»0 

DO  10U*1 ,h 

00  10K>1 ,N 
1*1*1 

lEU.EO.KlGj  >0  1? 

R I IL  l«0. 

GO  TQ  13 
H I (l  )•! . 

10  CONT INUE 
RETURN 
END 

Subroutine  NPR0(AtB,R.NRA,NCAR8,L.N.NCBI 
NPR3  PERFORNS  THE  MATRIX  PRODUCT  IAI*(B). 
RETURNS  RESULTING  MATRIX  INTO  (R>. 

called  from  umofx.pqpi 

DIMENSION  Am.BdURm 
LN»N  *L 

DO  1C  1*1, NRA 
DO  IC  J*1,NCB 
0*0. 

DO  10  K«  1 ,NC ARB 
IKA* (1-1  l•NCARB*K 
K J8»  |M*1  )*n:s*j 

1 JR*  tI-1 »»NCP*J 
0*0*  AI IK A)*a  CK JB  ) 

10  RllJR)>a 
RE  TURN 
END 


SUBROUTINE  C AR ( AR , L ABL ,NP  I 

CAR  PRINTS  OUT  THE  COMPLEX  MATRIX  (AR|  NAMED  (LABLl 
complex  ARIll 
NJ  ^NPPNP 
mR1TE(6,7)LAPL 
7 FORMATIIH  ,1A?| 

HR  IT  E I6.P  M AR ( J) . J*1 .N^  I 
6 F0RMAT(10E11.3> 

RETURN 

END 


SUBROUTINE  RVEC ( VEC * L ABL .NP ) 

RVEC  PRINTS  PUT  THE  VECTOR  (VECI  NAMED  (LABLl 

called  from  nline 

DIMENSION  VECCll 
HRITE(A,7)LABL 
7 FORMATIIH  ,1R2I 

HR  ITF(6,6MVEC(JI  .J*l  .NP) 

A FORRATIIH  >I0E11.3) 

return 

END 


IfSJ 


.1 


\ 


IHf  f IK0S((*S1 .311  ,H{CCS.HZSl««.(TC3S.f  rsIlk.Ml 
C THii  iUetauMNt  PtBfCtmi  l*«T{li'«fDlAft  CAIC  JlATiakS. 

C CAll(3  FADa  umOFi 

C JNPLF 1 HKCS.HZSlN.FTCOS.Frsll.&J.ijA.tJN.'iOA.FH.hZFIFt 

C JAA'J*t/KFL/»  ,IAP 

C3aA0*«/FPAB/A1,A^,A 

DATA 

0112 -CIL /2. 

sa«:MPt((o.,t.) 

lA •« /V 

Apsi  •psi/SFw^sre 
CPSI  <CS(AP&|  I 
Pa  *1 a*CP  si 

PP-(r»PK 

PN-IA-PA 

SP«SIII<PP«3I12I 

SH*S  ltllP«*0IL2l 

{ jK*C(>P(-Sj*XA*0XL2 l»SP/PP 

IF (pM.fe.o.i&oTa  11 

fjA'CfXPISJ^XAAaXl^l'SI/PA 

03  ro  12 

It  f J“-CtxP(SJ*n*3xi2l 
12  H3N>  H/f  AM  (a  ,A  ,AI  ,A2  • I AP  I 
Eti»i::**c*iPt--sj*P"«i;.xt2i 
H2cjs»r-«#{tjA.f  j«i 
a2SI  «»«-S  JA-F  JR  I 

FTC3S«M2C0S*2*CPS1 
f TSIS*M2SIA*2«CPSI 
«f  t JAN 
FN3 


SUSAOUTINF  AAA(AA.IABL.NP) 

C AAA  paints  out  THf  AATRti  (AR|  NANfO  (LAtLI 

C (AlLFU  FAQ*  POPl.NLlNl 

C I A| KSIJN  AA I I ) 

N2»NP*AP 
• AITFI6.71LAB1. 

7 FCA<AT(1H  .1A2> 

«A  IT  f (^,A  MAAIJ)  .J>1  .N2> 

0 FQAAATI ICFll .31 
AtTjAN 
FNO 

SjeAOUTINE  MAPlMT.E  .A.A1.A2.INPI 
COAAnM/rrOU/ei . B2.B3.BA.es. 66.B 7 
IF  ( IAPM.3.A 
3 E«l. 

AETjAN 

A TT«AL3C(A2/A1 >/lA2~Al ) 

TA-EXM  -Al*Tn-tXP|-A2*TTl 
AA  »A  /T*i 

E*AA  aIEXPI’AMT  )>EXP(-A2*T  n 
AETjRN 
2 CONTINUE 

eB*(  BS-96)*n  .•B3«T|*E|P(.B3*T> 

P9«CP6*B7)*U.*B2*TI*|xp(-e2*TI 

F10«-«S*(1.*B1*T1*EXP|-Bl*n 

Pll--e7*(l.*BA*T|*EXP(-es*TJ 

E*A«  (B0*B9*BIO«BI1  I 

AETjAN 

END 

SJBTnjTINE  CaaOO  I A ,9  .C ,Ni 

C CAA3D  A30S  lAi  N BT  N CJAPlFX  AATAICES  IAJ  AN3  (Bi. 
c called  faca  f in .UHOFX ,NL INE 

COMPLEX  Am ,E m .c< n 

00  2 J>1 .N 
2 (UI>AIJI*BlJl 
AETJAN 


Subroutine  saltia in.v .n.aout i 

C salt  AULTIPLIES  the  AATRIX  (A|NI  B»  TmE  CONSTANT  |VI. 

C CALLED  FAOA  F | N . P JP 1 ,NL I N E 
01  AE  NSION  RIM  1)  .ADUT  m 
NN«S*N 

CO  10  J*I.NN 
10  AOJT  IJI ■ V*A IN  I J) 

AETJAN 


SUBROUTINE  OELOtA.A.N.N.EPS. lEA) 

C &EL&  PFAFOANS  AATAlX  INVEASION  FOA  AEAL  AATAICES. 

C CALLED  FROM  P3PI 

C TAAENFRON  IBA  SYSTEA/3bO  SCIENTIFIC  SUBAOJTINE  PACKAGE.  1967. 
C REPORT  NO.  H20-01B6-<,.  PAGES  179-lBl. 

DIAE  NSION  ATI  )»Rm 

.»  DATA  2EAn/0.0002.HALF/O.S002.QNE/1.000/ 

|F<At23.23.l 
I IEA*0 

P1V2ER0 
AA'A*A 
NA-X  *A 
DO  3 L*l  .AA 
TB>A(L) 

IF (T P.LT  .2f AO)  TB*-TB 
IF ITB-PIVI3.3.2 


i-  24 


Ji 


i P J V T» 

l»l 

) CUNT  INut 

f 31  ‘ IPS***  I w 

DU  i T K> 1 ,■ 

U(Plvi 

« I F I Uk  u .*>«  7 

^ If  IP  IV-TIU  U.6,f 

k 118-K-l 

r p iv  I >u«it /«■  1 1 

j»u  -n/" 

I . I - J**-H 
J» J.  I -n 

00  J I -<»  .NH,» 

LI  *1 

TE  ‘P  !wl*»ILL  I 
aiLL  i 

• MU  • IP 

If  (4  . IB  t to 

If  Ij  I U .l/.io 

10  II*J*«* 

03  1 1 L>LSr.LI>lU 
TS>ML  i 
ll -I  •!  I 
*11  ; -*UL  > 

11  *llk  )»TB 

12  03  I » I >l  SI 
IL  *L  ♦! 

T«  «P  IV  l**UL  I 
AILL  )^*(L  I 

13  *IU  •TB 
*ILSn«J 

IST-IST*1 

J»0 

DO  16  II  >LST,UMC 
P U1  «-A(  in 
1ST* 

J«  J*  1 

Dj  IS  I « I ST 
LL»L-J 

All  ) »A(L  )«Pm  *4(11  1 

ra*A  (i. ) 

IMTP.lt  .ze  AO  I T6*-TE 
If  (TB-PlVllS,!^*!** 

lA  P(V>TB 

I *L 

IS  CCNTIfiUE 

OJ  16  L*K,NA,X 
LL*L  *J 

It  RILL  )=P(LL»*PnfI»R<L) 

17  LST»LST*M 

18  IFM-n  ?3,ZZ.|9 

19  IST«MNfH 
LST»M*l 

DO  21  1»2 ,H 
I [ »LST-1 
IST»IST-LST 
L«IST-P 
L'AimHALF 
00  21 
TB*<  m 
LL«J 

00  20  K*  I ST  ,*<M  ,H 
LL*L  L*l 

20  rd*TB>A(K)*R(LLl 

J*L 

RIJI  *R(K  I 

21  RIKI  *TB 

22  RETjRN 

23  IER>-I 
RETJRM 
END 

SU0R  OUT  1 NECOFLOIR , A ,M  *N,EPS, lER I 
C CGELG  PERFORMS  COMPLEX  MATRIX  INVERSION. 

C CALLED  FRQN  NLINE 

C 

c 

C3MPLE  X R ( I)  *A(  1 1 ,P  I VI »CTB 
IF (M  123, 23,1 
C 

C SEARCH  FOR  GREATEST  ELEMENT  IN  MATRIX  A 

1 1ER»0 
PIVO. 

MM  <M  *M 
NM -H  *M 
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•» 


00  ) L s 1 
T8=CAaS(AUI  I 
IM  T B-P  I VI  3,3  ,? 

2 P|V*Tb 
l*L 

3 CUNT  INUE 

C A(I)  IS  PIVOT  tUM^NT.  PIV  CONTAINS  THI  ABSOLUTE  VALUE  OF  A(1J. 

T3L=fPS®PlV 

C 

c 

C START  ELlHINAtlDN  LOOP 

L ST^  1 
DO  1 7 K=  1 
C 

C TEST  ON  SlNtULARITV 

IFIP  1V123,?3,<. 

<*  lFnf«)7,S,7 

5 I F IP  IV-TOL  Jb, 6, 7 

6 IERsK-1 

7 P I VI  =1  ./A(  I ) 
j»(  I -n /M 

J»J» 1-K 

C MR  IS  ROW-INDEX.  J«K  COLUNN-INOEX  OF  PIVOT  ELENENT 

C 

C PIVOT  ROW  SEDUCTION  AND  RQW  INTESCHANGE  IN  SIGHT  HAND  SIDE  S 
DO  3 L®KtNM,N 
LL=L  ♦! 

CTB=PIVI  ®R(LL  ) 

RILL  )=^R(L) 

8 RID  =CTB 
C 

C IS  ELININATION  TERMINATED 

1F(<(-M)9,18.16 
C 

t COLUMN  INTERCHANGE  IN  MATRIX  A 

9 LEND=LST*M-K 
1F(J)1?,1?,1D 

10 

DO  11  L*l ST .LEND 
CTBsAIL  > 

LL=L  ♦!  I 
AIL)  =A|LU 
11  MLLl^CTB 
C 

C ROW  INTERCHANGE  AND  PIVOT  ROW  REDUCTION  IN  MATRIX  A 

\2  DO  13  L =LSr ,MM,M 
LL^L  *I 

C Td=  PI  V I ®AILL  I 
AILL  )=A|L) 

13  AIL) =CTB 
C 

C SAVE  COLUMN  INTERCHANGE  INFORMATION 

XsREALIA(LST)  ) 

X*J 

AILSTlsCMPLXlX.O.) 

C 

C ELEMENT  REDUCTION  ANP  NEXT  PIVOT  SEARCH 

P I V=0. 

LSr=LST*l 

U*0 

00  16  I DLST  ,LEND 
P I VI  »“A I I I) 

IST«  tl«M 
J»J*  1 

00  IS  L > I ST  ,MM,M 
LL«L“J 

A ID  «A  IL  I ♦PI  V loAlLD 
TB«CABS(AID  I 
IFirB-PlVII5,l5,l4 

14  P I V»  T8 

1 *L 

15  CONT INUE 

00  16  L«K  ,NM,M 
LL*L+J 
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16  RUL  l«R(LLi«PlVl«R(U 

17  LST»LST*M 

C fMD  GF  FLtMlNATION  LOOP 

C 

C 

C PACK  SUBSTMUTION  AND  BACK  INTFKCHANCi 

18  lF('i-l)^3,^2,19 

19  lSTsNH*H 
LST»KM 

DO  2 1 1 a? ,H 
I I *L  ST-I 

isr=  tST-Lsr 

L'lST-- 

L = RE  AL ( A (L  n ♦ ,S 
00  2 1 Ja  n »NH ,M 
C TB=  RU» 

LL=J 

00  20  M *1  ST , 

LL*LL*1 

20  C TS=CTP-A  IK  J *R ILL ) 

K = J*  L 

R < J)  I 

21  « (K»  =CTB 

22  RETURN 
C 

C 

C ERROR  RETURN 

23  ItR»-l 
RETURN 
END 


1 


COMPLEX  FUNCTION  HZ F R E E ( W * A » A 1 , A2 » 1 M P I 
COMPLEX  B8,B9,UL0.Bll 

C HZFREEsCOMPLEX  AMPLITUDE  OF  INCIDENT  MAGNETIC  FREE  FIELD 

C IN  THE  TRANSVERSE  PLANE  OF  THE  NLINE  AS  A FUNCTION  OF  M(RAO/SEC). 

C CALLED  FROM  FIELDS 

C AsPEAK  TIME  DOMAIN  ELECTRIC  FREE  FIELD  AMPLITUDE  IN  VOLTS/METER 

C A1,A2=TWG  EXPONENTIAL  PULSE  SHAPE  PARAMETERS 

C IMP  SPECIFIES  TYPE  OF  MAGNETIC  FIELD  INPUT  PULSE  SHAPE 

C 1MP=0,  HZFREEaUNIT  IMPULSE 

C IMP= ♦! .HZFREE =TW0  EXPONENTIAL  PULSE 

C IMP* -I  .HZFREE  sSPEC  lAL  «»-E  XPONE  NT  I AL  PULSE 

C UMM  DN/YODLL/Bl , B2 ,B 3 , B9 ,BS , B6 ,B 7 
TO  =ALnG(A2/Al)/(A2-AlI 
TMaEXPt-Al*Tai-EXPI-A2*TOI 
XA=A/TM/377. 

IFUMP)S,2,3 

2 H2FREE=CHPLXn.,0.I 
GO  TO  9 

3 H2FREE=XA*n./CMPLX(ai,Ml-l./CMPLXEA2tWlJ 
GO  TO  <* 

5 XA=A /377 . 

P8s(B7*B6I*n.*02/CMPLX(B2.WlI/CMPLX(82,Ml 

B9a-BS»{l.*Bl/CMPLX(Bl,wn/CMPLX(Bl,M| 

PlOs(BS-B6)*(l.*03/CMPLX(B3,Kn/CMPLX(B3.M) 

B1 la -67* I I .♦P^/CMPLX(BA,M I I /CMPLX  tea  *« I 
H2FREFaXA®(8e*B9»BlO*Bin 
9 CDNT INUE 
RETURN 
END 

SUBROUTINE  HE  AO (1 » J » F R EK ,C XM G • I XAT ,V XMG . VX AT  »C OM NT ,T L ,P$ 1 , 1 Q • 
INM IRFS  ,I DNTFN.NP .NCOUNT,X) 

C MEAD  PRlNtS  OUT  THE  FREUUENCY  DOMAIN  RESULTS. 

C CALLED  FROM  NLINE 

0 I ME  NS  I ON  FREK 1600 i « C XMG I 6 00 f 1 0 I |CXAT(600»10I • VXMG ( 600 *I 0 I t 
IVX AT (600 ,101  ,C0MNT(8 ) 

C 

C 1*0,  PRINT  vx. .. 1*1  , PRINT  IX 

C J*PRINT  EVERY  J-TH  VALUE-U*!  PRINTS  AIL  FREK  VALJESl 

C 

tONTFN*IONTFN*l 
X X /TL 

WRITE  (6,600  II0NTFN,(C0MNT(RN),>(N>1,BI,IDNTFN 
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ll. 


Vx  ,wr ~.i  »v>i 


■t. 


600  F3R'4AT(1HUSH1.0.>*19,IOX8A10,10XSH1  .()•>,  19//) 

IMl 

3 URirF(6,b01)XX,U«PSI 

601  FURMATllH  , 1 <.  hVUL  T AOl  S AT  X * , IF  <»  . ? , 1 1 Nf  L t NO  TH  , 3 X/HL  * , 

lit  19  .9,1X<.HPSI*,1F6.2//) 

60A  FQRHATUH  •19HCURHFNTS  AT  X = , IF  9 • 2 • 1 ^H*L  I NE  L E NG  Th  » 3X2HL  « » 
lit  19  .9,3X9HPSI  = ,1F6.?//I 
hRITF(6,60S) 

605  F3R1ATUH  ,2%  /HRAD/SE  C * 1 X I iHCCiNOUC  TQR  1,TX|H?,10X1H3,10XIH9,10X1H5 
I .lOX iHh, lOXlH? .lOXlHP  ,10XIM9  ,9X?H10/  ) 

DD  10  K=1 ,NCOUNT ,J 

WR  IT  E(6 ,60?) FRFK (K  I , (VXNG(K,H ) ,H«1 »NP ) 

10  NRirE16,603)(VXAT(K,N),Hsl,NP) 

60?  FQR«^AT(IH  , 1 E 10  . ? , 1 1 E 1 1 . 3 ) 

603  FGRRATdH  , 6X  3H  AR  G , 1 X 1 1 1 U . 3 ) 

GO  TP  9 

? RR  IT  E16,609) XX, TL  ,PSI 
WR1TE(6,60S) 

00  11  Ks  1 ,NCQUNT  • J 

RRITE(6,6  0?)FREK(K), (CXHG(K,N)  ,H*I  ,NP ) 

11  WR1TF(6,603)(CXAT(K,N),H-1,NP) 

9  RETURN 

END 


SUBROUTINE  CSMLT(RIN,V,N,ROUT  ) 

C CALLED  FROM  F I N . UMQF X , NL I NE 

CO'iPLF  X ttlNU  ) ,ROUT  ( 1 ) ,V 

C MATRIX  QR  VECTOR  RINU)  MULTIPLIED  BT  A COMPLEX  CONSTANT,  V. 

00  10  J»l  ,N 
10  R3JT  (J}sV«RIN< J) 

RETURN 

END 


SUBROUTINE  FIN(X,BETA,PVX,PIX,Y,PI,NP) 

C This  SUBROUTINE  PERFORMS  INTERMEDIATE  CALCULATIONS. 

C called  FROM  NLINE 

DIME  NS  I ON  PI  I 100)  ,X  1 1100)  ,P8  (100  ),vnOO) 

COMPLEX  PVXU),Plxn),SJ,S,CC,XC(100),XXnOO),PPUOO),P9|lOO) 
SJ*CMPIX(0.,1.) 

N?*NP«NP 

CALL  XIDENT(NP,Xl ) 

DO  10  J«1  ,\2 

P9IJ  )*CMPLX(P1(J)  ,0.1 

10  XX (J  )^CMPLX( X 1 ( J ) ,0. ) 
the  T A=PE  TA*X 
C«C3S(THETA) 

S*SIN(ThETA)*SJ 

CALL  CSMLT(P9,S,N?,PVX  ) 

CC*CMPlX(C,0.) 

CALL  CSMLT(XX,CC,N?,XC) 

CALL  CMAOOIXC ,PVX ,PVX ,N?) 

CALL  SML  MPI  ,C  ,NP,P8  I 
CALL  C SMLT( XX  ,S,N?  ,XC 1 
03  1 1 J>1,N/ 

11  PP(J  )«CMPLX(P8( J)  ,0.) 

CALL  CMADD(PP,XC,XC,N2) 

CALL  5MLT(r,-l.,NP,xn 
DO  1 ? J*1  ,N2 

1?  XX  (J  )>CMPLX(X  n J) ,0. ) 

call  CHPRO(XX,XC,PIX,NP,NP,0,0,NP| 

RETURN 

END 

SUBROUTINE  SMAG(A,B,ACMT) 

C SNAG  CONVERTS  THE  COMPLEX  NUMBER  lA)  TO  POLAR  FORM.  B*MACNITUDE  OF  (A) 

c and  (acmt)  is  the  argument  OF  lAI. 

C CALL  ED  FROM  NLINE 

COMPLEX  A,C 
C>C3NJG(A) 

C»A*C 

? 0>SaRT(REAL(O) 

ARcREALlA) 

AI>A  IMAGIA) 

IF(AR)3,9,3 
9 IF(A  1)5, 6, 7 
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6 ACNTcO. 

cu  ro  6 

S A&m  *2  70  . 

GO  TO  ft 

7 AGNF  «90. 

GOTO  B 

3 ACHT *S7.29576*ATAN2(A1 (Aft) 

I F (AGHT  )9  ,S  ,B 
9 AGHT -AGftT*3bO • 

8 KETJRN 
END 

Soe^OaTlNE  Df  UTXU  ,N»4.ft| 

C CALLED  EAOA  POPl 

C DELETES  I-Th  COLUMN  AND  KDIi 

DIMENSION  AdUBd) 

N1*N>1 

L«0 

DO  10  J*lrN 
00  10  K*1.N 
IF  (J  .ED. I IGUTU  10 
IFU.EO.DGQTO  10 
L*L*  1 

ft(L)  *A(N*U-1  1«K1 
10  CONT  INUE 
RETJRN 
E NO 


OF  N*N  MATIII  A*  AHUKETUPNS  N1*N1  B. 


SUBROUTINE  UWUFX(X,tf,Z»PSI«UA.MX,V,N,ELj| 

CALL  FO  from  NLINE 

DIMENSION  Zd)tELJdO).CAPdOO)»ECM10),vm»EC2(10)«OdOO) 

CQMPLE  X UXd  dWXd  dSJtH7C0S,HZSIN,E  VCOStCE  L JUO  I .E  L UC  d 0 1 . 
lVYdr0),ECICd0).EC2Cd0).WBd0dE8(10l,EVSIN 
CDHMON/XSAV/CAP 
U*l. ?Sft637CbE -6 
SJ*CMPLX(0.,1.I 
N2*N*N 

CALL  FIELOS(PSl,X,HZCOS,HZSlN,tYCOS,ETSlN,#l 

HZCOS»MZCOS»SJ*M 

HZSI  N*HZSIN*li 

E»COS»EYCOS*SJ*M 

E VSl N*FYS1N*W 

DO  19  I « 1 «N 

C*0. 

03  9 J*1 ,N 
Na*N®(I-l)^J 
C*C*EL J( J)»CAP(NU) 

ECK  d*-C/U 
DO  10  J* 1 ,N 

ELJ:  (JdCMPLXlELJlj:  n.i 
CALL  CSMLT(ELJC,HZC0S,N,CELJ1 
00  1 3 J0*1  ,N2 
13  Q(JS)*Z( JOl 

CALL  MPR0(0,EC1 .eC2,N,N, 0.0,11 
DO  1 2 K*1 .N 

EC2C (K)=CMPLXCEC2CK1,0.1 
C<LL  CSMLTCEC2C ,E YSIN.N.EB I 
CALL  CMAODiCELJ.EB.UX.N) 

CALL  CSMLUELJC.HZSIK.N.ELUC) 

CALI  C SMI T IEC2C .EYCOS.N.EC IC I 

CALL  CMADD(ECIC.ELJC.WB,N1 

DO  11  J«1  ,N2 

YYI J l*CMPLX<YI Jl.0.1 

CALL  CMPROIVY. HB. NX, N.N. 0,0,1) 

return 

END 

SJBROUTI NE  P3PI  iNWiRES  ,IU,YO  ,Y1  ,KQ»PI ,V,Z,C) 

this  SUBROUTINE  PERFORMS  INTERMEDIATE  CALCULATIONS. 

called  from  NLINE 

dime  NS  I ON  POT  1)  ,P  I d 1 ,VOd  ),YI  (1  ),Y(  n ,z  m ,C  d I ,RI  UOO)  .CAP  dOO) 
1,CC( 100) .COPTIOO) 

COMMCN/XEL/Y,IHP 

common/xsav/cap 
data  EP/l.E-10/ 

NP=YWIRES‘1 

N2«NP*NP 

L»0 


10 


12 


11 


29 


ntk 


DU  iO  ,NwIkrS 
03  10  Js  I 1 Rt  S 
L*L*  1 

K0-NU1RES«(1-1I*J 
NR=SW|RES*(J-U*I 
CC  UOl^C  U) 

10  (C (SR1>CC(NQI 

C CC^EXPANOeO  CAPACITANCE  NATRIX 

CALL  RAR(CC«?HCC,NHIRES) 

CALL  XI0ENT(NP»RI) 

CALI  UTLMX(1U,NWIRES,CC.CAP) 

C CAPsREDUCEO  CAPACITANCE  NATRIX 

CALL  RARfCAP^^HC  tNP) 

CALL  SHLT(CAP,V,NP,n 
C Y»Y  matrix 

00  iO  KH»1  ,H2 

30  COP(KM)-CAP(KM| 

CALL  CrLG(RI.COP»NP,NP,LP,IER) 

C R1»CAP  INVERSE 

IM  1 ER  .NE  .OIGOTO  lb 
GO  in  17 

16  WR  II f(6,b01  II  .lER 

bOl  FORMATIIH  .AenUNSUCCE SSFOL  INVERSION  OF  CAPACITANCE  HATR IX*oN*  • 1 I 
13,3KAhIER**1I3////I 

17  C*l./V 

CALL  SMLTtRI,Q,NPWI 
C 2>Z  matrix 

CALL  RAR(RI,?HC-.NP1 
CALL  RAR(Y,2HY  .NP) 

CALL  RARIZ*2H2  ,NP) 

CALL  MPRDIZ»YI,P1 ,NP,NP,OtOtNP) 

CALL  MPRDU»Y0,P0tNP»NP,0»0»NP| 

RETURN 

END 


SUBROUTINE  HT TM < C ONNT , I ON T FN , VT , C T , T , NT  I ME , N I 
C HTIM  WRITES  TIME  DOMAIN  HEADINGS  IN  THE  PRJNT3UT. 

C CALLED  FROM  NLINE 

DIMENSION  VTm»CTm»T(l} 

DIMENSION  CQMNTCBI 

CDMiCN/XY2/<TCNT,NTH,TL,PSI*X 

COMMON/XFL/V,IMP 

tau*tl/v 

IDNTFNsIONTFN*! 

XX*X  /TL 

MR  IT  F(C  «600l ICNTFN.ICOMNT (MN I • MN« 1 »B 1 1 IDNTFN 
bOO  FORMAT  ( 1 HUSH  I .O.*.  19, 10X6A10,10XSH1  .O.s»19//) 
NRlTE(F,60nXX,TL,PSI 

601  FORMATIIH  ,1AHRESP0NSE  AT  X « , 1 F A . 2 , 1 *L I NE  L E NG TH • 3 X2HL « , 

1 IE  1A.A,3XAHPSI * ,1F6.2//I 
WRITE(6,606IV,TAU 

606  FORMATIIH  .30HCABLE  PROPAGATION  VELOCITY,  V * , 1 PI E9.2 , 1 HM£ TE RS/S E C 
1 ,2X35HTRANSI T TIME  OVER  ONE  CABLE  L E NOTH  * , 1 E 9. 2 , AMSE C .// I 
WR!TF(6,60SINTH 

605  FORMATIIH  , AX9M T I ME  I SE C U 7X5HV0L T S ,9XAHAMPS , 5X 
116HAT  CONDUCTOR  NO. ,12/1 
WR  IT  EI6,60A  1ITIJUVT(JUCTIJUJ*UNMME,N) 

60A  FORMATI 1P3E1A.A) 

return 

END 

subroutine  Y0T1IY0,YI,20,ZI,N) 

C CALL  EO  FROM  ML  INF 

OIMtNSIQN  VOID  ,Yim,£OI  1 1,2111  ),AI  1100  l,AQI  tool 
C SET  21  AND/OR  20  NEGATIVE  ILESS  THAN  2ER0I  TO  REPRESENT  OPEN 

C CIRCUIT. 

00  5 IR*1,N 

DO  5 UC*1R,N 

NQ«N*(  IR-n*JC 

NP«N*I  JC*n»IR 

IFI2  UNO)  .LT.O.IGOTQ  6 

At (N0)>1 ./2I INO) 

GD  TO  7 


30 


6 AUNQ|*0. 

7 If  (MMNQ  ) .L  r .0.)  GOTO  6 
AQ(N0)*l ./Z3(N0) 

GO  T n q 
e AO(NO)  >0  . 

9 An'iPl'AUHUI 
S AU(SP|cAD(NQi 
DO  10  I > 1 »N 

00  10  J* I ,N 

i-n*j 
NP«S*(J-1 »♦! 

Yl  *0  . 

V2»3  . 

1 F { I "J ) 1 , 2 • 1 
2 00  1 1 JJ«1*N 

JK>I 

H-N*  CjK-1  )*JJ 
Yl»r  UAl  (M) 

11  V2»Y2*A3(M1 
Y M'401*V  1 
Y0(N0I»Y2 
GO  TO  12 

1 Yl  (YOM-AUSai 
YU(MQ)>-AQ(NQ1 

12  Yl  (NP) »Y 1 (NOl 
10  Y0(NP)-YP(NU) 

«£  r JRN 
FND 


SUBROUTINE  TRANS(FW>FT,ri,TT, IU,IT,H,KDNTRL1  TRANS002 

C TRANS  PERFORNS  THE  NUHERICAL  FOURIER  TRANSFORM  OF  FREQUENCY  DOMAIN 

C TO  T IMF  DOMAIN. 

C CAUfO  FROM  NLINE 

C Fi^sARRAY  OF  Flhl  POINTS  COMPLEX  INPUT  TRANS003 

C FT*  ARRAY  OF  F|Tl  POINTS  OUTPUT  TRANSOOA 

C N»  ARRAY  OF  ri  POINTS  INPUT  TRANS005 

C TT*  ARRAY  OF  IIMt  POINTS  INPUT  TRAHSOOfr 

C U*  HUMBER  OF  N POINTS  TRAN5007 

C IT*  NUMBER  OF  TIME  POINTS  TRANSOOB 

C KDNTRL  * .NE.O  PRINT  BH * W , F ( U I . A ,6  ,C  • .EQ.O  NO  PRINT  TRANS009 

DIMENSION  FTIlO00)*TT(lOOO),H(lO0O)»A(SOOI,BCS0OI»C(SOO>  TRANSOIO 

DIMENSION  &(3)  TRANSOM 

COMPLEX  FWIIOOO)  TRANS012 

DOUBLE  PRECISION  T , C 1 ,C 2 • C 3 • C4 ,C S , C6 »C S2  .C SO  • SN2 »SNO • XP • X TRANS013 

DOUBLE  PRECISION  N 1 , M2  • , H 1 2,  H 22  , H 32 , ETl , E T 2 ,E  T 3,06  L T A TRANSOM 

GU)*2.  TRANSOIS 

G(2)*24.  TRANSOM 

G(31  *720  . TRANSOM 

C ARE  DIMENSIONS  EXCEEDED  TRANSOM 

IF  (I  H.GT  .1030  .OR.IT  .GT  .1300)  GQTOl  TRANSOM 

PI>3  .M1S927  TRANS020 

IFUONTRU  b.<»»6  TRANS021 

b WRITE(b»6001  TRANS022 

k I1=IW>3  TRANS023 

JJ»(  n/2  ) *2  TRANS02A 

IFIJJ.EO.in  1I*M«1  TRANS02S 

DOS  I*UIlf2  TRANS024 

Hl*4(n  TRANS027 

W2*M  (Mil  TRANS02B 

N3*W  (M2  ) TRANS029 

M12*  NM«2  TRANS030 

H22*H2««2  TRANS031 

H32*H3**2  TRANS032 

E TMRE  AL  (FW(  I M TRANS033 

E T2*  REAL ( FM( I ♦!})  TRANS0  3A 

ET3>REAL(FM(M2n  TRANS03S 

DEL  T A*  M?2*H 3 'H 32 »W2-W12*H3« W32*HMH  12 •W2-N22»H1  TRANS0  34 

L*I/2^1  TRAN5037 

A(U  *((M2-M3MtTl-(Nl-W3MET2^(Nl-N2MET3l20ELTA  TRANS0  3B 

B(LI *(-(H22-N32| *E T 1 ♦ ( W 1 2 -N32 » *6 T2 -( W1 2-N2 2 I *8 T3  1 /O E L TA  TRAWS0  3f 

C(LM(  (M229m3-H32»H2  MeTl-(ul2*tf3-W32*WlMET2«(Ul2*W2-W22qWl  MET3  TRANSOAO 
1 l/OfLTA  TRANSOAl 

IF  (K  0NTRU3,5 ,3  TRAN50A2 

3 BH*2 (I )»H/2.99E *8  TRANS043 

WRITECbrbOn  L • BH  , H (I  ) »F  .<  ( II  » A (L  > . B ( U ,C  ( L I TRANSD44 

601  FQRiATdS.FlO.b,  |P2E1S.3,E11  .3,3EIS.31  TRANS066 

600  F0RMATnHl,8X2MBH,MXlHN,13X2HRE  t9X2H|M,MXlMA,l6XlHB»MXlHC///l  TRANS066 

5 CONT  fNUE  TRANS067 

DOGD  ITC*1#IT  TRANS066 


I 


X-0. 0 

T-ir  n Tc ) 

037  1=^1, 

IF  ( ( h(  1)  *n  .LT  .1  .01  COTlXiI 
I F (i*  ( ) «T  ,CF  .3.it  7)  GaTU<»0 
C1»3(U/T**? 

.*A(l 
Ci*:  «Li/T 
C<»*^  .•AU)/I®*j 
Cb  = 3 K 1/ T 
Cb-A  (U/T 

cs^»ocrsun*^)*T) 
c sc=  ocns («( 1 1 «rj 
sn4'=dsi\(w( 

OSIN (ta( II »T  ) 

XP  s?./Pl*ICl»(CS<:TCS0l*C2*rt«(l*2)*CS2“«(lP 
l*CS»(»i<I*2|»SN2-W(l)*SN0)*C6*lMI*2l**2*SN2- 
0313  7 
kl  ril^D  (1  I 
H2  ( I *2  I 
AA=A(L) 

Bfi=3  (LI 

cc  ID 

Ql«ril*(AA»Wl**2/3.  * BB*Hl/2.  « CCl 


TRANS049 

TRANSOSO 

TRANS051 

TRANSQS2 

TRANS053 

TRANS0S4 

TRAMSOSS 

TRANSU^b 

TRANSOM? 

TRAfiSOSB 

TRANS0S9 

TRANS040 


CSO)  «(C3>C<>I*(SN2'SN0)TRANS06S 
Will *«2*SN0) ) TRANS066 

TRAKS067 
TRANS04B 
TRANS069 
TRANS070 
TRANS071 
TRAM507? 
TIAIIS073 


Q2*i<2«(  AA«ri2**2/3.  « BB*U2/2.  ♦ CC  ) 

SJN*0. 

If  It*  I®T  ,LT.  C.l  ) GO  TO  7i 
DO  7 1 K»U3 
XX  : < 

S=(-KI#*H*(».l*T|••(^•^)•wl•(AA*^^l••2/C3.♦^.•^XI♦BB•Wl/(2.♦^.•XX| 
X *CC/(l  .♦2.*XXI  1/GlKl 
71  SU“*SUK»S 
73  T£iiil*xuM*Jl 
SJ'^^0, 

I F ?®T  .IT.  O.l I GO  TU  03 
DO  b 1 K = 1 ,3 
XX  sK 

S = (-1.1**K*(»»2*T|*»(2»KI  •w2*  { AA*W2  •*2  / (3  .*2  . *X  X ) *88  ® h2/ ( 2 . +2  .*XX  | 
X *C  C/ ( I .*2 . *XX 1 1 /G(K I 
81  SUH*SUF'*S 
B3  TFK82=SuN*y2 

XP  = 2 ./PI » lTfc8N2-Tt  1 1 
7 X=X»XP 
<»0  F T ( I TC  I *X 
rftjrn 

1 >i8ITF(«>,bm  UfIT 

610  FCiRXAT(3<»H  OlMfNSlCNS  OF  TRANS  ARE  E XC  EE  DE  0 . //// <»H  1W*,I6,9H  IT*, 
I I6////UH  CALLING  OUHP./  I 
STOP 
END 


TRAN5076 
TRANS07S 
TRAN5076 
TRANS077 
TR  ANS076 
TkANS079 
TRANSOBO 
TRANSOBl 
TR  ANS082 
TRANSO03 
TRANS08<i 
TKANS08S 
TR  AN SO 86 
TRANS087 
TRANSOeS 
TRANS089 
TRANS090 
TRANS091 
TRANS092 
TRANS093 
TRANS09<» 
TRANS09S 
TRANS096 
TRANS097 
TRANS099 
TRANSIOO 


1 

i 

I 

I 


1 


7.  INPIT  DATA  KORMAI 

A sunimarv  ot  ihc  order  and  Idrmal  bv  which  data  arc  read  inlo  MdNL  is  ^iven  in  fORIRAN 
slaiements  below. 

READ(5.5000)A.A1 , A2 
READ(5 . 5001 )PERCEN 
READ(5.505)(COMNT(MN)  MN  1 . 8 ) 

READ(5.500)MSETS  NWIRES 
READ(5.501 )(CMAXWL(J) , J-1 , NCAPS ) 

READ(5 , 503) (ELI( J) , J=1  NP) 

READ(5 , 504 ) IQ 

D088  JJ88=1 .MSETS 

READ(5,501 )(ZI( J) . J=1  N2) 

READ(5 . 501 ) (Z0( J ) . J=1 ,N2 ) 

READ(5 , 502)PSI .TL 

88  CONTINUE 

5000  FORMAT (3E 10  0) 

5001  FORMAT(F5  0) 

505  FORMAT(8A10) 

500  FORMAT (215) 

501  FORMAT(8E10  0) 

504  FORMAT(I5) 

502  FORMAT(2F10  0) 

where.  NCAPS  = NWIRHS*(NW|Ri:S+  1)  2,  NP  = \MRHS-I.  and  N2  = NP*NP. 

A listing  of  a sample  data  set  consisting  of  iO  data  cards  follows.  The  FORTRAN  card  columns  are 
indicated. 

Data  Card  Columns 


5 

1 

10 

1 

20 

1 

30 

1 

40 

1 

50 

1 

60 

1 

1 

1 1 

1 

+ 2 

1 

3 +6 

1 

5 +6 

1 

1 

1 

1 

1 

1 

7 

1 

1 

1 

INE 

STANDARD 

DECK* • 

1 

1 

1 

1 

3 

1 

1 

1 

1 

+0 

1 +0 

1 +0 

3 279-  1 1 - 

1 704-11 

3 279-1  1 

3 74 

-7  1 

58-7 

1 . 

+2 

1 +2 

1 +2 

+2 

1 

+2 

1 ^-2 

1 +2 

+ 2 

90 

30 

8.  OLTPLT  FORMAT 


NLINE  output  in  the  Standard  Deck  is  in  line-printer  form.  The  user  may  add  optional  plotting 
^ instructions  at  the  locations  indicated  in  the  deck.  fTequency-domain  results  are  printed  out  when  the 

subroutine  HEAD  is  called.  Time-domain  results  are  printed  out  when  the  subroutine  HTI  M is  called.  The 
control  parameter  J in  the  listing  is  u.sed  to  limit  lines  printed  if  other  forms  of  output,  such  as  tapes  or  plots, 
are  used.  Informative  comments  are  included  in  the  listings  of  the  two  subroutines. 

■5  Figures  3 through  7 show  samples  of  a typical  NLINE  pnntout.  This  pnntout  is  generated  when  the 

NSD  IS  run  with  the  sample  data  listed  in  section  7. 

Figure  3 shows  the  first  page  of  the  output,  where  several  arrays  of  interest  are  printed  in  linear  form 
(i.e.,  the  ZI  array  shows  values  of  ^)-  ^re  the  terminating  impedance  arravs  defined 

in  section  4.  LE  is  the  magnetic-field  coupling  parameter  vector.  Yl  and  YO  are  defined  by  equations  (19) 
and  (20).  CC  is  the  Maxwell  capacitance  coefTicient  array.  C is  the  Maxwell  capacitance  array  after  the  / 
row  and  i columns  have  been  deleted  from  CC.  the  / ^ conductor  being  the  reference  conductor.  C"  is  the 


V.  33 

% 


Li. 


matrix  inverse  of  C.  Y and  Z are  defined  in  equations  (15).  PO  and  IM  are  defined  in  equations  (17)  and  (18), 
respectively. 

Figure  4 sliows  the  second  page  of  the  output  for  this  run.  fhe  radian  frequency,  amplitude,  and  pha.se 
of  the  voltage  induced  on  each  conductor  (see  listing  for  subroutine  UFAD  for  printout  instructions  of  the 
current,  if  desired).  F'or  e.xample.  the  first  frequency  is  1.1  X 10^  rad  s and  the  voltage  magnitude  on 
conductor  I is  4.'M  x 10  ^ V.  the  phase  is  given  directly  below,  marked  by  "ARC”,  and  has  the  value  of 
270  deg.  Similarly,  in  column  "2”,  the  voltage  and  phase  induced  on  conductor  2 is  given  as  2.09  x 10  at 
270  deg.  respectively. 

Figures  5 and  6 show  printouts  of  the  time-domain  voltages  and  currents  induced  on  conductors  I and 
2,  respectively.  The  heading  indicates  that  the  computations  were  made  at  .v  = 0,  on  a transmission  line  of 
length  / = 30  m,  with  angle  of  incidence  = 90  deg.  Cable  propagation  velocity  and  transit  time  over  one 
cable  length  are  also  pnnled. 

Figure  7 shows  printouts  of  the  incident  electric  field  used  for  the  computations.  The  exact  time-domain 
value  of  the  field  is  obtained  from  subroutine  TKMPLT,  and  is  printed  in  column  "EXACT.”  The  column 
"TRANSFORM”  gives  the  results  from  the  Nl.INf:  numerical  transform  of  the  frequency-domain  electric 
field  obtained  from  subroutine  HZFRF.F,.  fhe  close  agreement  for  most  of  the  values  indicates  the  relative 
accuracy  of  the  numerical  transformation. 


•*NL  lft(  STASDtRS  Df 


2 1 
20 

TO 

r} 

•20Df-Cl 

cc 

.103i*3i 

C 

C- 

t 

2 

.2fSF«0l 

.29Sl*31 


.icct«01 

.I33f 
.ISJf -Ob 
103E-OI 
-.lOCt-Ol 

-.I7D£-10 

.2UE*M 

.1C3EO) 

.nn-zi 


.I00f*03 
.loot *33 

-.lOCf-31 
-.lO&F-O I 
.1D0{*31 

.2  1 7E • I ; 

.103EO3 

.T63E-01 


.I00£*03 

.100E»33 


.?0P=-31 

.2C0C-01 

.1001*31  .328E-10  -.170f-l3 

.328e-|0 
.<•101  *1 1 
.b8«E-32 
*33 
*31 
.2<»se*31 


.IOCE*01  -.170E-10  .32«F-13 


Figure  3.  Sample  NI.INE  output  (first  page). 


34 


SttSOA&O 


VCLTiif^  4T  X»0.3J*U\t  LtHt'H  L« 


.iOOOt  *32  PS  I*  90.09 


nao/stc  CUHOuCn^  I 

.1  if  ,<.9«.c  - >9 

• 9&f*CS  .l6it*07 

.6U»9s  .nn-07 

AAO  .^63M0) 

.abt*9S  .391‘-J7 

*<C  ..“ftsMUJ 

.Uf»C.b  .S34t-07 

4«C  .7fc7{*J) 

.37f*JS  .lfe7t-CS 

i^C.  ,7t.U»3i 

.6?E*9b  .?77t-cb 

UC 

.«7E*0b  .3{mt-3t 

4Ai 

.I2f»37  .‘.liE-Vb 

AAC  .^«3E*03 

•37f*07  .ICftt-OS 

AAC  .ta7i*03 

.6^t•37  .lO^E-OS 

AA&  .IS3t*03 

.87t»37  .8371-56 

AA3  .131£*G3 

.UE*58  .6<.U-CJ6 

AAG  .IUE*C3 

.37E«0a  .1661-06 

AAC  .S3U*0; 

.62(*5S  .113E-06 

-AAG  .8S6C»02 

.87E*0S  .170f-C.7 

AA&  .3S^C*03 

.12E»0>  .SS6£-t7 

AA6  .SS6E*02 

.37E*39  .I48t~07 

AAC  .<i06E*32 

.62E*09  .399£->b9 

AAG  .S8^t*02 

.•7F«09  .29S£'0B 

A8&  .327C»03 


Figure  4.  Sample  NLINF:  frequency-domain  output  (second  page). 


••f.L  iNf  STANCAAO  5fC»-» 


2.10fc*C6Mt  I EA  S/ScC  t T<A\SIT  TI»>E  3VE8  C6E  (ABLE  LENGTH*  l.A3E-07StC. 


AT  C3'lDu:T3h  NO.  1 


S.7l',f.E-0  3 
A .7AtSE -02 
7.661  ‘!.02 
' .>A.t  • 

7.  7‘-«A-  -G. 

7.1  TiA{ 

6 .AS2AE -J? 
1 -02 
l.l7..-,f-';2 

U5>SSE  -3? 
7.Sl6AE-0» 
A.272CE-G3 
-I  . 7 76AE -02 
-I  .81366-02 
•1.71ASE-02 
-1  .6)301 -02 
-1.7065E-02 
-2.S3G5E-02 
-2.3635E-52 
-2.I387E-02 
-I.9309E-02 
-I  .9127E-02 
-2.1 I9AE-02 
-3.9370E-02 
-F.7169E-02 
•1 .S3So£-02 
-l.A9A3£-02 


Eigure  5.  Sample  NI.INH  time-domain  output  (third  page). 


1.0. • ^0ft400)  •••Hint  )T4ti0AR0  OtCK*^ 

RESPONiE  AT  l•C.OO•U'«E  LER&Tk  l>  .)000E*0? 
CARLE  PRORACATION  VELOCIIf.  V>  2 .1 OE •06Nt T (R S/ SEC . 


PSI*  90.03 

TRAdStT  TIME  OVER  OME  (ARlE  LEMCTh-  l.ASE'OTSEC. 


1.0. • A0A630) 


TintlSECl 

1 .3000E  >17 

2. A^8U-C9 

R. TIRIE  >C9 
2.271SE -OP 

S. W8TE  -OP 

T . 9RSAE -3A 
i.3®«.)E-07 
I.J700E-07 
1 .A*5  7E  -07 
1.9A1AE-07 
2.2272E-07 
2.9129E -07 
2.798AE  -07 
).0A^3E  -07 

3. J700E  -C7 
3.A9S7E-07 
i.iAUE  -07 
A.2272E-07 
A.S129:  -07 
A.790RE  -07 
5.0RA3E-O7 
S.3703E  -07 
S.AS57E-07 
5.9*l<.e-07 
6.2272E-07 
R.S129E-07 
A.7RR6E-07 
7.0RA1E-C7 


VOL  n 

-1.<.S3AE-0I 
-I.2092E *00 
-1 .9SJSE  *00 
-2.l23Pfc*00 
-1 .9«2M  *00 
-1.9A«.3E*00 
-l.A/8SE*03 
-l.b)2(;E*jO 
*3761-01 
-3.<.4A<.t-0l 
-3.ne«.E-ci 
-2.A029E-01 
-1 .4I27E-01 
A.J30SE-01 
L.2BieE-0| 
6.3SB6E-01 
3.R7SOE-01 
A.1I90E-01 
6.3728E-01 
S.9318E-01 
S.4170E-01 
L.AA39E-01 
6.883SE-01 
5.&2A5E-01 
R.0S33E-01 
6.A998E-01 
6.0223E-O1 
).«61AE-0i 


AMPS 

-6.7A7AE-0A 
-S.i9S0t-O3 
-9.3»6  ?l -03 
-9.9AS6E-C3 
-9.2b29E-0> 
-9.C79AE-03 
-6  .93061-03 
-7.63191-03 
-1  .leJiC -03 
-6.268SE-09 
-9  .4992E -C9 
-l.8  9<,lf-09 
1.980U-09 
2.6562E-03 
2.6993E-03 
2.98911-03 
2.3927E-03 
2 .3916E -03 
3.3910E-03 
2 .8691E-03 
2.9678E-03 
2 .3287E-03 
2 .2381E-03 
2.3090E-03 
•2.1099E-03 
1.8399E-03 
1 .69SAE-03 
I.SS9AE-03 
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Figure  6.  Sample  NLINE  lime-domam  output  (fourth  page). 
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6066003 


INCIDENT  ELECTRIC  EIELO  ERCN  EXACT  EjoATlOh  AND  FRO*  F CJA I E R TRAN9EJRN  IN  VUTS/NETER 


riNMSK) 

E aAC  T 

TAAN9FURN 

1.3C0CE  -12 

2 .6S96E -09 

7 .0/ JOE  *00 

2.898  IE -0« 

9.9667E*C1 

9 .710  IE  *01 

9.7I93E  -0« 

» . 1 IBbE  *01 

9 .1 PP  U *01 

2.27J9E -0* 

9. .6961 *01 

9 .9  3661  *0  1 

9.1287E  -OP 

9.9y66E*Ji 

9.9  76»F  *01 

7.i89?E -OP 

9.0830E  *01 

9.J9761 *01 

1 .0861E  -07 

8 .6b3U  «01 

8 .6X99E  *01 

1 .3700!  -07 

7.796U*CI 

T.7702t  *01 

1 .999  7E  -07 

7. I091E  *31 

7.3P6PE*0I 

1. 9616! -07 

6.6J23?«01 

6.607P1 *01 

2 .2272E-07 

9.7907f  .31 

9 .796IE  *01 

2.9129E-07 

9.I692E.0I 

9.1609F *31 

2.7986?  -07 

6 .99231 *31 

6 *31 

3.3P93E-07 

6 .3692F  *0 1 

6 .36091 *0  1 

3.3700E -07 

3.9P69E.0I 

3 .9602E  *01 

3.6997E  -C7 

3.1906E  *01 

3.1299F  *01 

3 .96161 -07 

2.T63')E«Ql 

2 .T  3601*01 

6.22  721  -07 

2.61221*01 

2.3877E*0I 

6.9129E  -07 

2.10261*01 

<.078  3E*01 

6.J9P6E  -07 

1 .82A7E  *01 

1 .80691  *01 

9.3861E  -07 

1.9d72f *0i 

1 .963IE  *01 

9.3700E -07 

1.97911*01 

1.39ltE*01 

9.6997E  -07 

1.1892E*0t 

1 .1  699E*0I 

9. 96161-07 

t.0269E*01 

1.0033E*01 

6.2272E  -07 

8 .A963E  *00 

8.6198E*00 

6.9129E-07 

7.6262E.00 

7.39121*00 

6.7986E-07 

6.9967E*00 

6.3299E«00 

7.3863E  -07 

9.6320E*00 

9.6C26E*00 

Figure  7.  Sample  NLINE  incident  EMP  output  (fifth  page). 
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9.  SAMPLE  PROBLEM 


This  sample  problem  is  compatible  with  the  \LINh  Standard  Deck  listing  and  with  the  sample  data  set 
included  in  this  manual.  I'he  procedure  lor  solving  a tspical  transmission-line  problem  with  NLINE  is 
outlined  in  this  section  and  computation  of  E:.MP  coupling  with  a configuration  of  three  parallel  conduits  is 
described.  The  F.MP  environment  and  the  environment  orientation  with  respect  to  the  conductors  are 
assumed  known. 

The  conductor  cross  section  for  the  sample  problem  is  shown  in  figure  8. 

For  the  sample  problem,  the  EMP  is  defined  to  arrive  at  broadside  incidence  90  deg).  Thus,  the 
magnetic  and  the  electric  fields  are  oriented  entirely  along  the  Z and  X directions.  TTie  Poynting  vector  P is 
onented  in  the  minus-y  direction,  as  shown  in  figure  8.  A specific  orientation  of  the  magnetic  field  with 
respect  to  the  conductor  cross  section  determines  a unique  set  of  values  for  the  magnetic-field  coupling 
parameters  (eq.  ( 10)). 

The  conductor  length  is  chosen  to  be  30  m and  all  terminations  are  chosen  to  be  100  ohms.  The 
propagation  velocity  among  the  conductors  is  chosen  to  be  70^?  of  the  free  space  velocity  of  light.  This  is  a 
typical  value  for  propagation  along  many  classes  of  transmission  lines  analyzed  in  nuclear  EMP  technology. 
The  point  of  computation  of  currents  and  voltages  is  chosen  at  the  ,v  = 0 end  of  the  conductors.  Pulse)  I)  is 
used  for  this  run.  The  thin-wire  approximations  are  u.sed  to  obtain  the  Maxwell’s  capacitance  coelTicients  and 
magnetic-field  coupling  parameters.^ 


E' 


Y 


COORDINATE  ORIENTATION 


Figure  8,  Sample  problem  conductor  cross  section,  and  onentation  with  respect  to  the  incident  EMP. 


2 Jams  Klehers,  Elecironwgnenc  Put^e  Cnuptm/;  wtih  Loxxte%s  XtulnainJut'ior  Transmission  lanes.  Harrs  Diamond  l.at>oratoncs  T R - 171 1 (June 
1975). 
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I'he  sample  problem  paramclers  are  summari/ed  below. 


(a)  IncidenlbAU’  l’iiKe(  li 

(b)  Angle  of  incidence:  V ‘*(1  den 

(c)  Onenlalion  of  IT/  as  shown  m figure  H 

(d)  Conduelor  length:  / = .'()  m 

(e)  Conductor  propagation  velocils : v - 2 1 x lO'’ m s. 

(f)  Conductor  terminations  (chosen  identicalK  at  ,r  = 0 and  .t  = /i: 
connected  as  shown  in  figure  I 

(g)  Reference  conductor  designated  as  the  /ero-th  conductor,  0, 

(h)  The  Ma.xw'ell's  capacitance  coelTicient  matn.x: 


100  ohms. 


1.*  l.*“ 

.12.8  — 1 7.0  pf  m 

-r.O  .12. 8_ 


•Arbitrary  capacitance  values  have  been  listed  for  the  reference  conductor,  since  NLINE  uses  only  the 
reduced  capacitance  matrix  for  these  computations: 


12.8  - 
-17,0 


The  full  capacitance  matrix  is  read  into  NI.INT.  allowing  the  u.ser  to  change  the  reference  conductor.  In 
that  case,  the  values  for  nil  capacitances  including  the  reference  conductor  must,  of  course,  be  included 
correctly  in  the  capacitance  matrix  (eq.  (25)). 

(i)  The  magnetic-field  coupling  parameters: 

I' 

“ “ Li  58J 

C ' and  T'  have  been  computed  bv  program  PTNTL,  listed  in  appendix  A. 

The  NLINE  Standard  Deck  accepts  all  the  above  parameters  and  performs  the  computations.  The  NSD 
printout  of  this  sample  problem  was  given  in  section  8.  figures  d through  II  show  optional  plotter  output  of 
the  results. 

Figure  9 shows  the  results  printed  under  “CONDLTTOR  I"  in  figure  4.  Results  in  figure  10  correspond 
with  the  "VOLTS"  column  in  figure  5.  Figure  11  shows  the  results  printed  under  "CON  Di  d OR  2"  in 
figure  4.  Results  in  figure  12  correspond  with  the  results  under  "VOLTS"  in  figure  6.  Finally,  figure  11  shows 
the  incident  electric  field  used  in  the  computation,  corresponding  to  the  column  labeled  "TR.ANSFORM"  in 
figure  7. 
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Figure  9.  SAMPLE  PROBLEM:  Frequenev -domain  plot  of' voltage  induced  on  conductor  No.  I by 

Pulset  I ). 
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riMf  <nO 

Figure  10.  SAMPLE  PROBLEM:  lime-domain  plot  of  voltage  induced  on  conductor  No.  1 by 

Pulset  1). 
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Figure  13.  INCIDENT  EMP  FOR  SAMPLE  PROBLEM:  Time-domain  plot  of  Pulse(  I)  obtained  by 
numerical  Fourier  transform  from  the  frequencv-domain  expression. 
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10.  NLINE  APPLICATIONS:  INTRODLCTION  OF  NONLINEAR  LOADS  BY  LLMPED- 

PARAMETER  NETWORK  SYNTHESIS. 

For  many  NLINE  computations  for  cables  in  Army  syster  s.  it  is  sufficient  to  represent  the  terminations 
of  the  multiconductor  transmission  line  by  “open."  "short,"  or  "matched"  resistive  loads  (i.e.,  load  resistance 
large  compared  to  the  characteristic  line  impedance,  load  resistance  small  compared  to  the  characteristic  line 
impedance,  or  load  resistance  close  to  the  characteristic  line  impedance).  Direct  computations  of  current  and 
voltage  are  made  by  NLINE  for  the  above  types  of  problems.  Nonlinear  loads,  such  as  terminal  protection 
devices,  surge  arrestors,  and  other  electrical  circuits  may  be  added  by  incorporating  the  NLINE  output  with 
circuit  analysis  codes— such  as  DAMTRAC— used  at  HDL.^  The  NLINE  code  represents  the  solution  for 
the  transmission-line  response  due  to  a coniinuouslv  distributed  EMP  excitation.  The  NLINE  output  is  not 
directly  compatible  with  typical  circuit  analysis  codes  in  this  form  since  the  circuit  analysis  codes  require 
discrete  lumped-parameter  circuits  and  discrete  current  or  voltage  sources  as  input  parameters. 

Figure  14  shows  symbolically  the  method  employed  to  convert  the  continuous  source  output  from  the 
transmission-line  theory  to  the  equivalent  discrete  voltage  source  and  lumped-parameter  network.  The 
approach  is  to  obtain  the  Thevenin  voltage  and  circuit  for  the  particular  transmission  line  and  for  the 
particular  EMP  illumination.  Once  this  has  been  done,  the  load  can  be  arbitrarily  specified  at  one  end  of  the 
transmission  line.  In  particular,  complex  circuits  can  be  attached  to  the  transmission  line  at  this  end  and  the 
response  can  be  computed  by  circuit  analysis  code.s.  By  this  method,  the  other  end  of  the  transmission  line 
must  still  have  real  loads  compatible  with  the  NLINE  Standard  Deck.  The  outline  below  summarizes  the 
NLINE  computations  required  for  the  network  synthesis. 


Baker.  A.  McNuU.  B Shea.  [).  Rubenstem.  Danuj^v  4nalvsi.s  '^lodtfied  1 R4C  Computer  Program  iDiMl RACt.  Harr\  Diamond 
Laboralones  TVi-75-6(May  1975) 
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EQUIVALENT  CIRCUIT 

nline  synthesis 


Figure  14.  Lumped-parameter  network  synthesis  by  NLINE  computations. 

NLINE  THEVENTN  EQUIVALENT  CTRCLTT  METHODS  FOR  USE  WITH  CIRCUIT  ANALYSIS 
CODES  AND  NONLINEAR  LOADS  AT  ONE  END  OF  THE  TRANS.MLSSION  LINE 

(a)  Specify  all  NLINE  input  parameters  for  the  particular  problem. 

(b)  Designate  one  end  of  the  transmission  line  (i.e..  .e  = 0 or  .r  = /)  where  complex  load  conditions  are  to 
be  analyzed. 

(c)  Specify  loads  at  the  other  end  of  the  transmission  line  which  best  approximate  the  actual  terminating 
conditions  and  are  compatible  with  the  NSD  (/.e..  real  loads  having  values  best  approximated  to  be 
large,  intermediate,  equal,  or  small  compared  to  the  transmission-lins  impedance). 

(d)  Leaving  the  above  parameters  fixed,  use  NLINE  to  compute  the  open  circuit  voltage  and  the  short 
circuit  current  at  the  circuit  end.  To  accomplish  this,  run  NLINE  once  with  arbitrary  negative  values  of 
load  impedances  at  the  circuit  end  (NLINE  treats  negative  terminating  resistance  values  as  open 
circuits).  Store  the  open  circuit  voltages  U.^iia)  on  tape.  Next,  compute  the  short  circuit  current  /^(u'). 
with  small  (e.;^.,  I ohm  or  less)  values  for  the  terminations  at  the  circuit  end,  and  store  on  tape. 


(e)  Compute  the  transmission-line  impedance  /^tw)  b\ 


/.ua)  = 


Construct  an  equivalent  circuit  with  lumped  parameters  for  /j(ui).  For  broadside  h.MP  incidence. 
1^  = 90  deg.  I I clearly  exhibits  the  impedance  characteristics  of  a senes  of  parallel  R-L-C  circuits 
whose  component  values  are  determined  from  the  Qualils  Factor.  Q.  which  is  obtained  from  ' Z.ue)  |.  Q 
is  determined  from  the  fundamental  resonant  frequency/,,  of  i /,(u;)  | and  the  corresponding  half-ptrwer 
frequencies  /,  and  /j  from*" 


Knowing  Q.  the  values  of  R.  /..  and  C are  determined  from 


C>  = R C 


The  equivalent  lumped-parameter  circuit  is  constructed  from  a series  of  R-L-C  circuits,  as  shown  in 
figure  14.  where  the  value  of  1.  is  varied,  holding  R and  C constant.  Generally,  a low-frequencv  element 
must  be  included  to  account  for  the  low-frequency  charactenstics  of  i /'a(w')  j.  For  the  following  example, 
this  component  is  a resistance  R^.  The  values  of  1.  are  chosen  so  that  each  lank  circuit  resonates  at  the 
odd  harmonics  of  /„.  Hence,  the  values  of  are 

L, 

Li  = .m,  = 1.3 (xJd 

OT,  ' (31) 

and  ( = 1.2 /VV.,  where  A',,  is  the  number  of  tank  circuits  u.sed  in  the  svnihesis.  is  determined  from 

Q.  R.  and  in  the  second  of  equations  (30).  The  larger  the  value  of"  .V,  used,  the  higher  the  high- 

frequency  cutofT  point  beyond  which  the  transmission  line  is  no  longer  adequately  represented  bv  the 
lumped-parameter  model. 

The  input  impedance  of  the  equivalent  lumped-parameter  circuit  approximates  very  closeK 

the  iransmission-lme  impedance  /„(u!)  for  frequencies  less  than  the  high  frequency  cutoff.  Hence. 
can  be  used  to  replace  the  continuous-source  transmission-line  impedance  by  a specified  senes  of  R-L-C 
tank  circuits  programable  into  circuit  analysis  codes. 

(g)  The  load  voltage  at  the  circuit  end.  due  to  an  arbitrary  load.  computed  from 

,,  f oc(w')Zl(u>) 

‘ (32) 

or  from: 


= Z',(u,)  -b  /,>■)  ■ 


fcquations  (32)  and  (33)  describe  Thevenin  equivalent  circuits  for  the  transmission  line 


(h)  A;,  a final  step,  l is  transformed  to  the  time  domain.  is  programmed  in  the  circuit  analysis  code 

along  with  the  lumped-parameter  representation  ol  /'„  and  the  time  domain  representation  of  I he 
problem  is  solved  in  the  time  domain  by  the  circuit  analysis  code. 

Ihe  Nl.lNF  Standard  Deck  was  modified  to  perform  the  above  steps.  The  sample  problem  data 
parameters  were  used  for  the  computations.  Figures  15  and  16  show  the  open  circuit  voltage  and  the 
short  circuit  current,  respectively.  The  transmission-hne  impedance  computed  by  NLINE  is  shown  in 
figure  17.  The  corresponding  result  for  the  lumped-parameter  model  appears  in  figure  18.  Twenty  R-L-C 
circuits  were  used,  giving  a high-frequency  cutolV  at  70  MFlz.  A higher  cutoff  can  be  attained  by 
increasing  the  number  of  circuits  used  m the  model. 

Equation  (33)  is  used  to  test  the  accuracy  of  the  model.  Figure  19  shows  the  voltage  on  Conductor  No.  I 
of  the  .sample  problem  computed  from  the  lumped-parameter  model.  Compan.son  of  this  result  with  the  result 
from  figure  10  shows  that  the  lumped-parameter  model  gives  a relatively  accurate  representation  of  the 
transmission-line  response.  The  frequency  and  time  arrays  are  selected  m the  NSD  on  the  basis  of 
transmission-hne  length  and  on  five  transit  times  of  the  pulse  as  it  is  reflected  back  and  forth  along  the 
transmission  line.  Some  of  the  high-frequency  spikes  appearing  in  the  lumped-parameter  results  are  due  to 
the  numerical  Fourier  transform  "noise,"  and  can  be  corrected  by  a choice  of  more  appropriate  frequency 
and  time  points  in  the  transform.  In  general,  an  optimum  frequency  and  time-array  selection  scheme  is  not 
available  for  a fixed  number  of  points  (600  in  the  NSD)  and  for  a wide  range  of  iransmission-line  lengths. 
Guidance  on  numerical  transforms  is  given  in  the  DNA  EMP  Flandbook,^  and  should  be  used  to  select  the 
frequency  and  time  arrays  appropriate  for  a .specific  problem.  When  properly  used,  the  lumped-parameter 
technique  provides  a useful  application  of  the  NLINE  code  for  problems  where  circuit  analysis  codes  must 
be  used  to  treat  E.MP  coupling  with  transmission  lines  terminated  into  complex  circuits. 

32 
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Figure  15.  Open-circuit  voltage  computed  by  NLINE  for  use  with  the  lumped-parameter  network  model. 


I,  Huichinii.  [iffense  Nuclear  Agcncv  fMP  HamlKnik  Chapter  14  PMP  Data  Analysis.  BDM  A- I?5-7S1  R-R 1.  24  Sept  1975 
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Figure  16.  Short-circuit  current  computed  by  NLINE  for  use  with  the  lumped-parameter  network 

model, 
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Figure  17.  Transmission- line  impedance  computed  by  NLINE—/  = 30  m.  i/- = 90  deg.  Q — 1.16. 
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Figure  18.  Pqurvaienr  lumpcd-parameier  network  transmission-line  impedance  obtained  from  20  R-I.-C 
circuits  / = 30  m.  = 90  deg.  Q — 1.19. 
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Figure  19.  Time-domain  plot  of  voltage  induced  on  conductor  No.  1;  computed  from  the  lumped- 

parameter  network  model. 
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II.  APP1.K  ATION  OF  M.INF  FO  HFAI.  ( A8I.F.S:  < OMPAHI.SO.N  OF  M .MKRK  Al. 

PRF.DKTIONS  WI  IH  MF.A.Sl  RF.l)  RFSl MS. 

The  NI.INE:  code  is  designed  to  predict  nuclear  or  sinuilator  generated  hMP  coupling  with  cables  used 
in  systems  under  FMP  evaluation  A conlldencc  level  must  be  established  for  the  N1  INF  computations  in 
order  for  the  code  to  be  successlulK  applied.  Ihe  .AF.SOP  radiating-dipole  LMP  facility  at  HOI.  was  used  to 
obtain  e.xperimental  coupling  data  for  comparison  with  Nl.lNh  predictions.  The  results  of  that  test  are 
discussed  in  this  section  following  a summary  of  important  aspects  to  be  considered  when  applying  NT.INh 
to  real  cables. 

The  NLINE  code  requires  several  input  parameters  discussed  in  previous  sections  which  either  must  be 
determined  experimentally  or  computed  by  approximate  methods  for  complex  cables.  How  accurately  these 
parameters  are  specified  for  the  experimental  conditions  ultimately  determines  how  well  Nl.INF. 
computations  agree  with  measured  results.  The  accuracy  in  measuring  dilTerential  voltages  or  currents  in  the 

radiated  EMP  experimental  environment  also  can  be  expected  to  vary  with  the  type  of  cable  tested  and  j 

instrumentation  problems  (c.g..  common-mode  rejection  capability,  low  signal  level,  etc.)  encountered  in  the  j 

experiment.  The  above  factors  can  be  expected  to  produce  ditferences  between  the  computed  and  measured  1 

results. 

Coupling  measurements  and  NLINfi  computations  are  presented  in  this  section  for  .1-m-long 
transmission-line  samples  made  up  of  three  and  eleven  conductors  illuminated  bv  the  AFiSOP  environment. 

Determination  of  the  coupling  parameters  for  these  cables  has  been  discussed.-  The  cables  w'ere  oriented  in  ] 

the  test  area  parallel  to  the  AESOP  antenna,  and  were  located  035  m above  ground  at  a range  of  75  m on 

the  center-line.  The  significant  differential-mode  coupling  is  due  only  to  the  radial  component  of  the  AESOP  j 

magnetic  field  in  this  configuration.  Figure  20  shows  the  radial  magnetic-field  component,  along  with  the  j 

Pulse!  I)  approximation  used  in  NLINE  for  the  computations.  The  general  waveshape  and  peak  amplitude  ! 

for  Pulse!  I)  adequately  represent  the  AESOP  pulse.  Fhe  early  time  ramp  and  the  irregular  plateau 

characteristics  have  not  been  included  in  the  curve  fit.  Some  error  results  from  these  omissions.  To  eliminate  ' 

this  type  of  error,  the  user  can  modify  COMPl.EiX  FUNCTION  HZFREE  in  NLINE  to  return  the  exact  ! 

frequency  charactefistics  of  the  measured  environment. 

I ! 

Unshielded  transmisstorr  lines  made  up  of  two  and  three  conductors  represent  an  important  class  of  I 

interest  for  nuclear  EMP  analysis.  A very  large  group  of  Army  systems  uses  these  types  of  transmission  lines 
for  communication  and  power-transmission.  The  NLINE  code  can  be  used  with  a high  confidence  level  for 
' these  types  of  cables.  Figure  21  shows  the  measured  and  computed  results  for  the  three-conductor  test 

sample.  The  peak  values  agree  very  closely,  but  more  ringing  appears  in  the  measurement  than  in  the  ] 

computation  of  the  load  current.  This  is  probablv  due  to  the  reactive  load  instrumentation  conditions  I 

existing  in  the  expenment  which  are  not  accounted  for  in  the  lossless  NLINE  model.  ! 

i 

The  increase  in  the  number  of  conductors  m the  transmission  line  increases  the  margin  of  expected 
difference  between  predictive  computations  and  measurements.  The  decrease  in  accuracy  for  determining  the 
coupling  parameter  values  and  increased  lossy  conditions  in  the  cable  can  account  for  this  difference. 

Figure  22!a)  shows  the  eleven-conductor  sample  cross-section  with  the  computed  and  measured  coupling 
^ parameters  listed  by  each  conductor.  Agreement  is  not  always  giKtd  between  the  experimental  and  numerical 

methods  employed  here.  Admittedly,  more  accurate  means  are  available  for  the  user  to  determine  the 
parameters.  The  thin-wire  approximations  are  used  here  to  demonstrate  the  degree  of  success  attainable  from 
/ the  easy  to  u.se  approximations.  Also,  for  some  cables,  values  obtained  by  these  methods  may  be  as  gotxl  as 

, those  obtained  by  more  exact  methods  in  the  practical  case,  where  the  coupling  parameters  are  subject  to 

J change  !due  to  cable  stretching,  twisting,  bending,  etc.)  with  every  deployment  of  the  cable. 

Figure  22!b)  lists  the  measured  and  computed  peak  currents  for  the  eleven-conductor  sample.  The  i//i 
i variation  is  computed  from  i 

2 Jams  Klebcrs.  fUectronuigncnc  Puhv  C nupUn^  with  /.dv5/«'\5  \fu/iu  unJut  tor  7 rttn\mission  l.invs.  Harr\  Diamond  I Ahi'raU'nes  1 R - PI  1 (June 
1975) 
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CURRKM  (nr 


L?:  THIN-WIRE 
' APPROXIMATION 

IL®)  : MEASURED  VALUE 


X 10®  H 


Cjj:  COMPUTED  FROM  POTENTIAL 
COEFFICIENTS 
(Cjj):  MEASURED  VALUES 


(a)  C i>nduLi«T  ^comctrs  and  cahlc-coupling  paramclcrs 


.Measured 

Computed  by 

'onductor  \o. 

Peak  Current  m m.A 

NLINE 

dB  Variation 

1 

+500 

+340 

+ 3.35 

2 

-400 

-475 

- 1.49 

3 

-240 

-205 

+ 1.37 

4 

-140 

-190 

- 2.65 

5 

+280 

+ 85* 

+10.35 

6 

+280 

+230 

+ 1.71 

7 

+340 

+330 

+ 0.26 

8 

+240 

+100 

+ 7.60 

9 

+120 

- 45* 

* 

10 

+120 

- 48* 

« 

|H|  FV.ik  mdiikcd 

».urrcntN  McaMired  and  n'nipuicd  for 

the  M S( >P  ensironmeni 

• Indicates  large  error 

intrtKiuced  eilher  in  the  measurement 

or  computation  of  the 

cou- 

pling  parameters. 


Figure  22.  Comparisons  of  measured  and  computed  results  for  the  eleven-conductor  transmission  line. 
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iHi  = 20  log, measured  peak  computed  peak)  , 


(34) 


It  IS  evident  from  the  table  (fig.  22b)  that  the  eondmons  for  conductors  5.  9,  and  10  were  not  adequatelv 
treated,  either  bv  the  experiment,  or  b\  the  computation.  However,  at  least  six  conductors  show  relativeh 
good  agreement  between  the  computed  and  measured  results. 

In  conclusion,  the  application  of  the  NI.INh  code  to  real  cables  requires  u.se  of  g(K>d  engineering 
judgement  in  order  to  obtain  reliable  results.  Whenever  possible,  it  is  desirable  to  combine  both  computations 
and  radiated  EMP  environment  response  measurements  for  the  cable  under  analysis,  \ level  of  confidence 
can  then  be  established  for  the  NI.INE:  computations.  It  is  then  possible  to  introduce  circuit  analysis  codes, 
to  change  terminations  or  other  NLINE  input  parameters,  and  use  the  NI.INE  code  as  a powerful  tool  in 
EMP  vulnerability  assessment  and  hardening  programs. 
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APPKNDIX  A.  PK(M;RAM  1 ISHN(;  KOK  I UK  ( APA(  I I AN(  K AM)  magnk  i k -hei.d 
K)l  Pl.INt;  PARAMKI  KR  ( ODK.  K/  V/7, 


3 

The  determination  of  the  capacitive  and  the  induciiNe  coupling  parameters  lt>r  mulliconduclor 
transmission  lines  has  been  defined  and  documented  bv  many  investigators.  st>me  of  uhich  are  listed.' 

The  methods  for  determining  the  parameters  vars  widely  with  respect  to  the  degree  of  accuracy  and  the  ; 

degree  of  ditTiculls  in  computation.  Generally,  all  the  melht>ds  are  assumed  to  have  the  same  degree  of 
accuracy  when  the  conductor  spacings  are  large  relative  to  the  conductor  radii,  and  the  thin-wire  conditions 

are  approached.  The  thin-wire  approximations  appear  to  be  adequate  lor  some  applications,  and  are  desirable  : 

because  of  ease  in  computing  the  parameters. 


The  program  in  this  section  computes  the  .Maxwelfs  capacitance  coefficients  from  the  potential 
coefficients,  and  computes  the  magnetic-field  coupling  parameters  from  the  ihin-wire  approximation.  The 
number  of  conductors,  .V.  and  the  reference  conductor.  IQ.  are  designated  by  the  user  in  the  main  deck  as 
indicated  by  the  *‘CU"  comments.  An  80-characier  FORTRAN  comment  COMNT  is  read  in  by 

R£AD(5 , 500) (COMNT(MN) ,MN=1 .8) 

500  FORMAT(8A10) 

The  cable-conductor  cross-section  coordinates  are  read  into  the  program  by 

READ(5 , 501 ) (R(NN) ,T(NN) , NN=1 .N) 

501  FORMAT(8F10  3) 

where  R(NN),T(NN)  are  the  r.O  polar  coordinates  of  the  centers  of  each  conductor  of  the  cable.  The  variable 
r is  in  meters  and  9 is  in  degrees.  The  angle  9 is  measured  clockwise  fiom  the  y-axis,  as  shown  in  figure  8. 
body  of  this  report. 


The  conductor  radii  A(l)  are  given  in  meters  and  are  specified  by  the  user  in  the  main  deck.  .As  an 
example,  the  coupling  parameters  for  the  Sample  Problem  were  computed  by  PTNTI..  The  listing  follows. 


PTNTI (INPUT .OUTPUT ,TAPf 5* INPUT ,T*PE6»0UT PUT  1 j 

CU  subroutines  CELG.  XICENT.  and  OELETX  from  NIUE  BUST  ht  i 

C AlTIC'lFD  T3  HIS  OEC’^.  .| 

c i 

C THIS  PRDGRA'^  COMPUTES  ThF  REOUCET  “4X»tELl"$  CAPACITANCE  HATRIX,  AND  THE  I 

C HAGNETIC  FIELD  COUPLING  PARAMETERS  FCR  THE  MULT  I CDNOUCTDR  TRANSMISSION  j 

C LINE.  The  CAPACITANCE  MATRIX  IS  COMPUTED  3Y  THE  POTENTIAL  CPEFFICUNT  ! 

C method,  and  the  MAGNETIC  FIELD  COUPLING  PARAMETERS  ARE  CCjMPjTEO  'RCM  THE  j 

c thin  mue  apprdx I ma r IuN . i 

DIMENSION  R1  (100)  ,PI  (Un  .PP  (U1  ) i 

DIMENSION  p(ii.ii).x(in,Yni).R(in,Tm).3(u.m.A(n),  | 

ICQMNT(A) 

data  C, rad/1. 8E10..C17<.S3*?R/ 

CU  USER  SETS  tPF=  TO  EFFECTIVE  VALUE  DE  THE  RELATIVE  DIELECTRIC 

c constant  of  the  caple. 

f PR»  1. 

CU  USER  SETS  N*  TO  NUMBER  OF  CONDUCTORS  IN  THE  CA9LE- 
N«3 

CU  USER  SETS  I3»  TO  THE  REFERENCE  fPNDJCTOR  NUMBER. 

IG‘l 

KREF  I 

CP  USER  GENERATED  COMMENT  OF  FORTRAN  CHARACTERS  READ  IN  FROM  COLUMNS  1 

c Through  po  on  data  card. 

REAO(S.500)ICOMNT(MN).MN=1.5)  I 

SCO  FORMATIBAlOl 

CP  COORDINATES  OF  THE  CABLE  CONDUCTTR  CROSS  SECTION  ARE  READ 
C here  , IN  POLAR  FORM. 

C THETA  IS  MEASURED  CLOCKWISE  FROM  THE  V-AIIS. 

C R IS  IN  METERS*  AND  THETA  IS  IN  DEGREES. 


* J D.  Clements  and  C R Paul.  ( ompuuiiionof  fhc  < apiuiiann'  for  l)u‘U‘<  iru  -(  iHiieti  H irv\.  I m\ervit\  ol  Ken  tusks . R \ IX  I R-'U 

89.  Phase  Report  (March  1974) 

*S.  Frankel.  Cable  and  Muhuondudor  Tranxnusxion  line  Analius.  Harrs  Dianu»nd  I aK»ra  tones  I R-09|.I  (28  June  I9'4i 
^ K Kajfe/.  Uuliunnductor  Tran.xmi^si«n  Interaction  Nt>ie  Air  Force  Weapons  I aboralors.  Kirtland  \l  B.  (June  !9"2) 

C'  Kevser.  Modeling  Techniques  for  Mullu'onduclor  i ahles.  Theors  and  Prattue  Air  T »»rce  Wcap«’ns  I aK)ralor\  AI  W I TR-'’2-89 
(March  1973). 

^S.  Ramo  and  J R W'hinners.  helds  and  W«m’\  m Minlcrn  Radio.  John  W ilcs  and  Sons.  Ini.  . 2nd  ed  . New  N otk  ( IN(4i  2Nl 
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ftMO<S,SOlMS  (NN  ) ,T(KN),NN>1  ,N)  I 

SOI  fOtiATISFlO.))  J 

f P«1  .I-IS  I 

DA*»0.  ',| 

CU  USM  PROViDtS  VALUES  FCR  AdJ'RAOluS  OF  ThE  I-TH  CONDuCTPR,  ,! 

C FDR  ALL  CCN0UCT3RS  IN  TmE  TRANSHISSICN  LINE.  '' 

A(U  *.U8l 

AE?I*.13S  - 

A I 3 J * . I D S 
C 3 1 1 J » 1 , N 

X < J)  sR ( J ) ®SI N (RADOT I J I ) 

11Y(J)»R<JI®C3S(«A0®T(J)I 

03  I ? 1 s I 

0 3 I ? J»1  »N  1 

IFU  .E0.J)G3TQ13  ! 

en»j)>scRr((x(i)'X(j)i**2«^(V(i)>Y(ji)**?) 

C3T312 

13  H( I*  J) *A { J)  ‘ 

I?  C3NT  INUE  ! 

03  lA  I»1  ,N  j 

03  I A J.l.N  j 

lA  P II , JI  *C  *AL3G  (B  ( I ,KkE  F ) *?  ( KRE  F , J )/6  ( ! , J I /B  ( KRE  F ,KRE  F M/EPR  J 

hR  ir  EC  A ,599MCDNNT  fNN)  ,E  I 

KRirf(6,70?) 

70?  FORRATIlw  ,9X1HX,9X1HY/) 

OC  1 6 Js I 

lb  MR  IT  El A.6Q0) ( X( j|  ,Y(J}  ) | 

MR  I T E (fc , 703  ) ] 

703  FDRRAT  nHO#??HPOTENT  lAL  COEFFICIENTS///!  1 

Ob  1 S I > 1 ,N  ! 

15  mRITECA,601)(P{1 ,M) »M=l,*n  i 

599  FORiATIlMUBAlO///) 

, 600  FDRRATIIH  ,2E1A.A>  ’! 

601  FCRRATIIH  .IPUE12.3J 
CALL  XrOENT(NP,RI) 

j l-O 

00  1 9 J* 1 ,N 

DO  1 9 K«1  ,N  1 

L-L*l  i 

19  PI  (.  I=P(  J,M  j 

CALL  OFLFTXI IO,N,PI ,PP 1 1 

’ CALL  GFLC IRI ,PP,NP,NP,EP , lER J I 

C R I »C  APAC 1 TANC  E NA  TR I X i 

L*0  5 

Du  1 7 J*1  ,NP  I 

00  17  K*|  ,NP 
L*L*1 

17  PCj.KJsRMU 

MR  IT  Elb  »70n  IER*DAA 

701  FQRS  AT  HH-,  AH  lERs  , 1 1 3,3X  ,1F6.2,5H»A(J»  ///)  j 

MR  IT  FI  A , 700  I ' 

730  F3RRAT(1H-,1PhCAPACITANCE  MATRIX/)  -j 

03  I?  1*1, NP  i 

18  MR  IT  E 16 ,601  M PCI  ,M  ) ,M>  1 ) 

MRlTF(6,e03» 

800  F3RRAT(IH0,7RHMAGNETIC  FIELD  COUPLING  PARAMETERS  L (11 , L ( 2 ) • . . . , L ( N 

11  , WITH  REFERENCE  L (101 *0// 1 ' 

ft  00  SCI  K 8*1,  N 

XL*1.256637E-6»(Y|Kb)-Y|IQ)J  j 

601  MRirEI6,802)XL  j 

802  A0RMATIIP1E17.7)  ! 

STJP  I 

END  ;j 
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ATTN  CJIAPXLS  vr. 

Am:  PICHAPJj  F.  7CH/\EPEP 

att;;  v;illI/j;  '.  KArE/xS 
Am;  PPARD  r.  schle  il 

RCA  CORJCPATION 

x?'/ew?;e!;t  & U':;umercial  eyftfms 

ASTRO  ELECTRONICS  CIVISI.'N 
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SCIENCE  AFFLICATIONS,  INC. 
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